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SUMMARY 


An  oil  spill  behavior  model  has  been  developed  which  describes  the  changing 
area,  thickness  and  physical  properties  of  spills  of  Prudhoe  Bay  Crude  Oil, 
No.  2  fuel  (home  heating)  o*'  and  No.  4  (heavy)  fuel  oil  under  Arctic  marine 
conditions.  The  oil  properties  calculated  are  density,  viscosity,  pour  point, 
aqueous  solubility,  flash  point,  fire  point  and  interfacial  tension.  The 
model  calculations  Include  spreading  Into  thin  (sheen)  and  thick  slicks, 
drift,  evaporation,  dispersion,  and  formation  of  water-in-oll  emulsions 
(mousse  formation).  The  model  was  calibrated  by  fitting  the  model  equations 
to  data  obtained  from  outdoor  weathering  experiments  conducted  at  the  Coast 
Guard  RAD  Center,  Groton,  Connecticut,  during  the  Winters  of  1979/80  and 
1980/81  in  which  arctic  spring  conditions  were  simulated.  A  comprehensive 
sampling  program  yielded  data  on  the  changing  oil  properties.  The  model 
successfully  described  the  rate  of  evaporation  and  the  physical  properties 
under  these  test  conditions,  with  the  exception  that  interfacial  tensions  were 
not  well  predicted.  This  exception  is  probably  due  to  the  formation  of 
surface  active  compounds  by  oil  oxidation,  and  water-in-oil  emulsion  formation 
induced  by  rainfall.  It  is  believed  that  the  model  provides  the  capability  of 
predicting,  with  acceptable  accuracy,  the  behavior  and  properties  of  the  three 
oils  when  spilled  under  arctic  marine  conditions. 
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1 .0  INTRODUCTION 


There  is  an  obvious  incentive  to  develop  a  capability  to  predict 
the  extent  of  weathering  of  an  oil  spill  under  Arctic  conditions,  and 
the  resulting  changes  in  the  oil's  physical  and  chemical  properties. 

These  properties  influence  both  oil  behavior,  such  as  spreading,  and  the 
applicability  and  effectiveness  of  the  various  available  countermeasures. 

The  ultimate  objective  is  to  be  able  to  respond  rapidly  in  the  event  of  an 
Arctic  oil  spill.  A  computer  program  which  calculates  oil  property  changes 
as  a  result  of  receiving  input  data  on  the  oil's  properties  and  amount,  and 
the  environmental  conditions,  would  assist  this  objective. 

Many  oil  spill  models  have  been  developed  to  assist  with  countermeasures, 
most  being  of  the  trajectory  type  in  which  the  aim  is  to  calculate  where  the 
oil  may  be  at  any  time.  Such  models  essentially  add  wind  and  current  vectors 
for  selected  time  steps.  In  the  model  described  here,  the  aim  is  to  predict 
property  changes  rather  than  location,  thus  the  model  is  complementary  to  a 
trajectory  model.  It  may  be  wholly  or  partly  incorporated  into  a  trajectory 
model  if  desired. 

It  is  also  obvious  that  some  set  of  oil  property  measurements  should  be 
made  to  permit  rates  of  processes  such  as  evaporation  to  be  calculated.  The 
algorithms  or  equations  for  such  calculations  are  not  always  available  and  of 
a  form  compatible  with  the  data  acquired  or  with  the  oil  spill  model.  These 
algorithms  usually  arise  from  controlled  laboratory  or  outdoor  experiments  in 
which  the  rate  of  a  process  is  studied  in  isolation. 

In  this  work,  we  have  addressed  the  problem  of  integrating  information 
from  (i)  physical  and  chemical  property  measurements 
(ii)  outdoor  evaporation  experiments 

(iii)  recently  developed  algorithms  describing  oil  spill  processes 
into  a  comprehensive  oil  spill  behavior  model. 

Components  (i)  and  (iii)  relied  heavily  on  previous  work  undertaken  at 
the  University  of  Toronto  with  Environment  Canada  (Arctic  Marine  Oil  spill 
Program)  support  as  described  in  several  contract  reports,  notably  Mackay  and 
Paterson  (1980),  Mackay  (1978)  and  Mackay,  Buist,  Mascarenhas,  and  Paterson 
(1979).  Measurement  and  correlation  of  physical  properties  were  made  at  the 
University  of  Toronto.  Component  (ii)  was  undertaken  at  the  U.S.  Coast  Guard 
RAD  Center,  Groton,  CT.  A  final  integrating  effort  involved  joint  analysis  of 
the  results  and  formulation  of  the  oil  spill  behavior  model. 

Three  oils  were  tested  both  in  the  laboratory  and  at  Groton;  home  heating 
or  No.  2  fuel  oil,  heavy  or  No.  4  fuel  oil,  and  Prudhoe  Bay  crude  oil.  These 
oils  were  selected  by  the  U.S.  Coast  Guard  as  being  representative  of  oil 
types  most  likely  to  be  encountered  in  spill  situations  in  offshore  Alaska. 

The  approach  taken  in  this  report  is  to  first  describe  in  Section  2.0 
the  oil  spill  model  (a  complete  listing  being  given  In  Appendix  A).  The 
experimental  weathering  work  and  the  physical /chemical  determinations  are 
described  in  Section  3.0  and  the  results  in  Section  4.0.  The  calibration 
and  a  discussion  of  the  results  are  given  in  Section  5.0. 
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2.0  OIL  SPILL  MODEL 


Oil  spill  models  are  of  two  general  types;  trajectory  models  In  which 
the  aim  Is  to  predict  the  location  and  size  of  a  slick  In  a  given  geographical 
setting  under  the  Influence  of  prevailing  wind  and  current  vectors,  and 
behavior  models  In  which  the  properties  of  the  oil  are  predicted  In  some 
detail  from  a  knowledge  of  the  oil  composition  and  environmental  conditions, 
especially  temperature  and  wind  speed.  These  models  may  become  computation¬ 
ally  complex  and  rarely  are  they  combined.  A  recent  comprehensive  review  of 
the  status  of  models  has  been  prepared  by  Huang  and  Monastero  (1980).  Several 
models  are  described  In  the  proceedings  edited  by  Mackay  and  Paterson  (1980). 
Valuable  reviews  of  models  and  their  component  processes  have  been  provided  by 
Jordan  and  Payne  (1980). 

Possibly  the  most  advanced  models  currently  available  or  under  develoment 
are  those  developed  for  the  U.S.  National  Oceanic  and  Atmospheric  Administra¬ 
tion  Including  recent  work  on  their  behalf  by  Science  Applications,  Inc., 
Klrsteln  and  Payne  (1982). 

To  be  useful,  a  model  must  be  applicable  to  different  oils,  preferably 
by  inserting  a  number  of  parameters  into  the  program  which  characterize  such 
things  as  evaporation  and  changes  in  properties  such  as  density  as  a  function 
of  extent  of  weathering  and  temperature.  The  model  equations  should  minimize 
risk  of  computational  failure.  The  number  of  parameters  should  be  minimized 
and  the  inherent  assumptions  should  be  apparent  to  the  user.  These  require¬ 
ments  favor  simplicity  but  not  at  the  cost  of  loss  of  realism.  The  model 
described  here  represents  an  attempt  to  achieve  this  goal  by  modifying  and 
updating  a  model  developed  for  the  Environment  Canada  Arctic  Marine  Oil  spill 
Program  (AMOP)  In  1975-1980  and  most  recently  described  by  Mackay  and  Paterson 
(1980). 

A  complete  printout  of  the  model  Is  given  In  Appendix  A  with  behavior 
and  property  parameters  defined  Individually  for  Prudhoe  Bay  crude  oil,  home 
heating  oil  and  heavy  fuel  oil.  The  principal  features  of  the  model  are 
described  below. 

The  program  Is  written  In  FORTRAN.  It  commences  by  setting  spill 
conditions  of  oil  volume  and  the  duration  of  the  release.  Provision  for 
Instantaneous  or  continuous  releases  is  available.  Environmental  conditions 
of  temperature,  wind  speed  and  sea  state  are  defined.  The  oil  properties  are 
entered  in  the  form  of  appropriate  constants  for  the  equations  tabulated  later. 

Rate  constants  for  dispersion,  emulsification,  spreading,  drift  and 
evaporation  are  calculated  or  entered.  Provision  is  made  to  allow  any  process 
to  be  stopped  In  order  to  explore  Its  effect  by  comparing  "with"  and  "without" 
outputs. 

The  Initial  oil  properties  and  spill  thickness  are  set  and  the  calcula¬ 
tion  proceeds  In  time  Increments  by  computing  the  amount  of  oil  evaporated 
and  dispersed,  water  uptake  which  forms  water-ln-oil  emulsions,  spreading  and 
drift.  A  complete  printout  then  follows  giving  the  mass  balance,  the  rates  of 
the  processes  and  the  newly  computed  oil  properties  which  have  changed  as  a 
result  of  evaporation.  Any  new  oil  is  then  added  (If  the  spill  Is  continuous) 
and  the  calculation  repeated. 


The  oil  is  assumed  to  spread  into  a  thin  sheen  surrounding  thicker 
patches.  Most  of  the  spill  area  (e.g.,  70  to  90%)  is  sheen,  but  most  of 
the  oil  volume  is  contained  in  the  thick  patches.  Equations  are  given  to 
calculate  the  increase  of  the  thick  and  sheen  areas  during  a  time  increment. 
Provision  is  made  for  the  thick  slick  to  "feed"  the  thin  slick,  a  process 
which  is  observed  on  the  ocean  surface.  The  spreading  equations  have  been 
described  by  Mackay,  Buist,  Mascarenhas  and  Paterson  (1979). 

The  evaporation  calculation  uses  the  “evaporative  exposure"  concept 
in  which  the  oil's  vapor  pressure  is  first  calculated  as  a  function  of 
temperature  and  fraction  evaporated,  then  the  amount  evaporated  calculated 
from  a  conventional  mass  transfer  flux  equation.  Separate  calculations  are 
included  for  sheen  and  thick  slicks.  The  physical  basis  for  these  equations 
has  been  described  by  Mackay,  Paterson  and  Nadeau  (1980)  and  by  Mackay  and 
Stiver  (In  preparation). 

Dispersion  rates  for  thick  and  sheen  slicks  are  calculated  as  a  function 
of  windspeed  and  oil-water  Interfacial  tension  and  are  a  simplification  of 
equations  described  earlier  by  Mackay  and  Paterson  (1980). 

Emulsion  (mousse)  formation  is  included  using  the  equations  outlined  by 
Mackay  and  Zagorskl  (1982).  The  water  content  and  viscosity  of  the  oil  are 
calculated. 

Finally,  the  location  and  shape  of  the  slick  are  calculated  using  a  drift 
factor  of  3.5%  of  the  wind  speed.  It  is  possible  to  overlay  the  slick 
location  and  size  In  a  hydrographic  chart. 

The  variables  used  In  the  program  are  defined  and  units  given  in  most 
cases  by  comment  cards  in  the  program  but  additional  information  is  provided 
in  Appendix  A. 

A  specimen  output  of  the  state  of  spill  of  Prudhoe  Bay  crude  oil  is  given 
in  Table  2.1. 
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Table  2.1  Specimen  Output  for  State  of  Spill 
Prudhoe  Bay  Crude 
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3.0  EXPERIMENTAL  PROCEDURES 
3.1  Weathering  Experiment 

A  laboratory-scale  outdoor  experiment  was  conducted  at  the  Coast 
Guard  RAD  Center  during  the  Winters  of  1979/80  and  1980/81  to  simulate  the 
weathering  of  oil  under  arctic  conditions.  Four  oils  were  weathered  outdoors 
on  water  In  open  boxes  under  natural  environmental  conditions  for  periods  of 
up  to  a  month  with  periodic  sampling.  Temperature  and  wind  speed  were  moni¬ 
tored  continuously  to  allow  calculation  of  hourly  "evaporative  exposure" 
values.  Solar  radiation  was  also  monitored  to  elucidate  its  effect  on 
increasing  the  oil  temperature. 

The  primary  objective  of  the  experiment  was  to  document  the  effect 
of  evaporation  on  oil  weathering  under  cold  conditions,  and  its  role  in 
determining  the  oil  fraction  remaining,  and  the  physical  and  chemical  pro¬ 
perties  of  the  oil.  Of  particular  Interest  to  the  Coast  Guard  with  regard  to 
the  physical  properties  were  the  density,  viscosity,  and  combustibility  of  the 
oil,  as  these  properties  are  likely  to  dictate  the  effectiveness  of  various 
cleanup  measures  such  as  skimming,  use  of  suction  hoses  and  burning. 

A  second  objective  was  to  determine  the  dependence  of  the  rate  of 
evaporation  on  the  "evaporative  exposure"  and  to  verify  the  accuracy  of  the 
empirical  formulations  for  modelling  this  dependence. 

A  third  objective  of  the  experiment  was  to  investigate  the 
dependence  of  the  oil/air  interface  temperature  on  the  air  temperature, 
underlying  surface  temperature,  and  solar  radiation.  Of  particular  interest 
was  the  role  of  the  solar  radiation  in  raising  the  temperature  of  darker  oils 
above  ambient  environmental  levels. 

This  report  deals  with  the  first  and  second  objectives  as  described 
above.  The  results  stemming  from  the  third  objective  are  discussed  in  a 
separate  Coast  Guard  RAD  Center  Report  (Tebeau,  Meehan,  and  Myers,  1982). 

The  oil  weathering  experiments  were  carried  out  on  the  roof  of  the 
main  RAD  Center  building  at  A very  Point,  Groton,  CT.  This  location  offered  a 
marine  environment  with  ample  exposure  to  the  elements.  It  was  hoped  that 
during  the  winter  months,  this  environment  would  approximate  conditions  on  the 
North  Slope  of  Alaska  during  late  spring  and  summer,  when  large  quantities  of 
oil  would  most  likely  reach  the  ice  surface. 

The  experimental  apparatus  consisted  of  plywood  weathering  boxes  (60 
cm  x  60  cm  x  30  cm  deep)  lined  with  5  cm  of  styrofoam  insulation  and  fitted 
with  a  polyethylene  liner  to  make  the  boxes  reusable.  Transparent  covers  were 
also  fabricated  and  mounted  35  cm  above  the  oil  so  that  the  weathering  oil 
samples  could  be  protected  from  rain  and  snow  if  desired  as  shown  in  Figures 
3.1  and  3.2. 

The  experimental  procedure  called  for  freezing  a  block  of  fresh 
water  ice  in  each  box  in  the  RAD  Center  cold  room.  The  boxes  were  then 
transported  to  the  roof  where  a  400mL  sample  of  the  oil  was  deposited  on  the 
ice  surface.  The  No.  2  and  No.  4  oils  were  obtained  from  local  University  of 
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Connecticut  physical  plant  stocks  at  Avery  Point.  The  Prudhoe  Bay  crude  oil 
was  obtained  by  courtesy  of  the  Atlantic  Richfield  Company.  In  warmer 
weather,  the  oil  was  initially  deposited  in  a  (30  cm  x  30  cm)  cavity  chipped 
in  the  center  of  the  ice  block  to  prevent  it  from  running  down  the  sides  of 
the  block.  Usually  after  a  few  days,  enough  melt  water  was  present  so  that 
the  oil  covered  most  of  the  surface  area  inside  the  box.  The  oil  samples  were 
then  allowed  to  weather  for  periods  ranging  from  approximately  1  day  to  1 
month.  At  specified  times  during  the  period,  a  sample  of  each  type  of  oil  was 
removed  from  the  weathering  boxes.  Descriptions  of  the  environmental  monitor¬ 
ing  instrumentation  are  described  by  Tebeau,  Meehan,  and  Myers  (1982). 

Altogether,  109  weathered  oil  samples  (30  No.  2  oil,  31  No.  4  oil, 

35  Prudhoe  Bay  Crude,  and  13  No.  6  oil)  were  collected  and  analyzed  for  the 
Winter  1979/80  and  Winter  1980/81  experiments.  In  addition,  unweathered 
standard  samples  of  each  type  of  oil  were  taken  at  intervals  throughout  the 
experiment  period.  The  time  weathering  began,  time  weathering  ended,  duration 
of  weathering  in  hours,  average  temperature  during  the  weathering  period, 
average  wind  speed  during  the  weathering  period,  and  cumulative  evaporative 
exposure  for  each  weathered  sample  were  calculated,  the  results  being  tabu¬ 
lated  in  Section  4.1. 

A  preliminary  physical  properties  analysis  was  conducted  at  the 
Coast  Guard  RAD  Center  to  check  mass  fraction  remaining,  specific  gravity, 
viscosity,  water  content,  flash  point,  fire  point,  and  combustibility. 

Chemical  analysis  at  the  RAD  Center  involved  running  Gas  Chromatograph  (GC) 
scans  on  each  sample  to  delineate  its  chemical  composition.  A  complete 
description  of  this  analysis  and  the  results  are  given  in  a  companion  report 
by  Tebeau,  Meehan,  and  Myers  (1982).  Based  on  the  initial  results  of  the  RAD 
Center  analysis,  various  samples  of  the  No.  2  oil.  No.  4  oil,  and  Prudhoe  Bay 
crude  oil  were  selected  for  further  analysis  by  the  University  of  Toronto  as 
described  in  Section  2.0.  It  is  the  result  of  this  more  detailed  analysis 
that  forms  the  basis  of  the  model  calibration  effort  described  in  Section  5.0. 

3.2  Physical  Property  Determination 

Detailed  analyses  were  conducted  on  13  samples  of  Prudhoe 
Bay  crude  oil  (designated  by  a  prefix  "PBC"),  13  samples  of  home  heating  oil 
(HH),  and  13  samples  of  heavy  fuel  oil  (F0),  most  obtained  from  the  outdoor 
weathering  experiment  at  Groton.  The  remainder  were  samples  of  the  fresh  oils 
and  oils  obtained  by  artificial  weathering  in  the  laboratory  (designated  AB). 
The  artificially  weathered  samples  were  prepared  by  bubbling  a  stream  of  air 
of  known  flow  rate  through  a  volume  of  oil  for  specified  times  generally 
following  the  procedure  described  by  Mackey  and  Stiver  (1982).  This  latter 
group  of  samples  gave  results  which  were  used  to  develop  equations  describing 
the  physical  properties  as  a  function  of  temperature  and  extent  of  evapora¬ 
tion.  The  equations  were  then  tested  on  the  Groton  samples. 


Density 

Density  was  measured  at  0°C  and  20°C  using  a  standard  pycnometer 
calibrated  with  distilled  water.  The  results  are  presented  in  units  of 
kg/m3.  Specific  gravity  is  thus  this  number  divided  by  1000. 


Pour  Point 


Pour  points  were  measured  using  the  standard  American  Society  for  Testing 
Materlals-Institute  of  Petroleum  Method  #D97.  The  general  procedure  was  to 
heat  the  oil  sample  to  approximately  46<>c,  then  cool  the  sample  steadily 
with  cold  water,  then  Ice  water,  then  an  acetone-dry  ice  mixture  (as 
necessary)  until  a  temperature  is  reached  at  which  the  oil  in  the  test  vessel 
does  not  flow  when  tilted  horizontally.  The  results  are  expressed  in  °C. 

Viscosity 

Viscosities  were  measured  using  a  standard  Canon  Fenske  viscometer  at 
0°C,  10°C,  20°C,  the  results  being  reported  in  units  of  centi poise 
(cp).  The  now-standard  unit  is  Pascal-seconds  (Pas),  1  Pas  corresponding  to 
1000  cp  or  1  milli -Pas  (mPas)  to  1  cp.  In  most  cases,  it  was  either  not 
necessary  or  it  was  impossible  to  measure  viscosities  at  all  the  above 
temperatures .  For  example,  it  is  not  feasible  to  measure  viscosities  below 
the  pour  point,  thus  the  Prudhoe  Bay  samples  which  had  pour  points  in  the 
range  -2  to  +18°C  were  examined  at  10  and  20°C  or  20  and  30°C.  Only  two 
determinations  are  necessary  to  obtain  information  about  the  temperature 
dependence  of  viscosity. 

Aqueous  Solubility 

Approximately  5-10  mL  of  the  oil  samples  were  added  to  a  separatory  fun¬ 
nel  containing  50-100  ml  of  double  distilled  water.  The  separatory  funnels 
were  shaken  gently  in  a  wrist  action  shaker  for  about  two  hours  and  then 
allowed  to  settle  for  at  least  48  hours  before  analysis  at  20°c.  2-3  mL  of 
the  saturated  aqueous  solution  were  drawn  to  be  analyzed  by  "purge-and-trap" 
gas  chromatography.  The  results  are  given  in  g/m3  which  is  numerically 
equal  to  mg/L  or  ppm. 

The  gas  chromatograph  was  a  Hewlett-Packard  Model  5840  equipped  with  a 
flame  ionization  detector  and  a  Model  7675A  purge-and-trap  sampler.  The  GC 
column  was  a  wide  bore  50  meter  long  glass  capillary  column  coated  with  SE-30 
(WC0T  from  Alltech  Associates,  Inc.). 

Typical  chromatograms  of  the  fresh  oils  and  field  samples  are  shown 
in  Section  4.0.  It  was  noted  that  the  dissolved  components  were  mostly 
aromatic  in  nature,  e.g.,  benzene,  toluene,  xylenes,  substituted  benzenes 
and  naphthalenes  which  is  expected  as  their  solubilities  in  water  are  higher 
than  those  of  alkanes.  The  solubilities  were  determined  using  benzene  for 
calibration  purposes.  As  the  oil  weathered,  its  solubility  decreased  leaving 
the  higher  molecular  weight  components  such  as  polynuclear  aromatics  in  the 
solution.  Since  the  solution  was  purged  at  room  temperature,  It  was  Impos¬ 
sible  to  strip  hydrocarbons  heavier  than  n-C]6  or  naphthalene  from  the  water. 

This  is  not  regarded  as  a  serious  problem  because  these  hydrocarbons  have 
low  solubilities  and  they  will  contribute  little  in  addition  to  the  measured 
solubility. 
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It  should  be  noted  that  these  solubilities  are  measured  In  distilled 
water,  not  seawater.  It  Is,  however,  relatively  simple  to  calculate  the 
seawater  solubility  from  that  of  distilled  water,  the  seawater  value  generally 
being  75%  of  the  distilled  water  value. 

Oil -Water  Interfacial  Tension 

Oil -water  interfacial  tension  was  measured  by  the  spinning  drop  technique 
at  20°C  in  which  the  shape  of  an  oil  drop  Is  observed  while  spinning  In  a 
horizontal  column  of  synthetic  seawater  (35  g/L  sodium  chloride),  centrifugal 
effects  causing  the  oil  to  locate  along  the  axis.  The  oil  drop  adopts  a 
nearly  cylindrical  shape  with  the  relative  length  and  diameter  controlled  by 
the  balance  between  centrifugal  forces  (which  tend  to  elongate  the  drop)  and 
Interfacial  tension  forces  (which  tend  to  favor  area  reduction  and  thus  resist 
elongation).  By  measurement  of  droplet  volume  and  length  at  known  spinning 
frequency,  it  is  possible  to  calculate  the  interfacial  tension,  if  phase 
density  information  is  known.  The  technique  has  been  described  by  Mackay 
and  Hossaln  (1982)  who  give  the  equations  which  are  used  in  this  work,  namely, 

(1/Y)1/3  (3X/2) (H)l/3  =  H  r3/Y  +  1 
where  H  =  Apui2 

and  y  is  interfacial  tension  (dynes/cm),  X  is  semi -major  drop  axis  (cm), 
i.e.,  half  the  length,  r  is  the  radius  of  the  drop  sphere  of  equivalent  volume 
(cm),  a)  is  angular  velocity  (s"')»  and  Ap  is  density  difference  between  oil 
and  water  (g/cm3),  The  equation  is  solved  iteratively  for  Y 

The  calculation  if  done  In  SI  units  gives  interfacial  tension  1q 
Newtons/meter  (N/m).  The  commonly  used  CGS  unit  of  1  dyne/cm  is  10"  3  N/m  or 
1  milli -Newtons/meter  (mN/m).  An  advantage  of  the  technique  Is  that  there  is 
no  solid  surface  in  contact  with  the  oil  as  occurs  in  the  pendant,  sessile 
drop,  or  de  Nouy  ring  tests. 

For  selected  samples  of  low  viscosity,  oil-air  interfacial  tensions  were 
also  measured  at  20°C  (using  the  de  Nouy  ring  method). 

Oil  sheen  spreading  properties  are  influenced  by  these  interfacial 
tensions,  the  net  spreading  tendency  (ST)  being  given  by 

ST  =  Y  AW  "  Y  OW  "  Y  OA 

where  the  Interfacial  tensions  (y)  are  subscripted  AW  (air/water  Interface), 

OW  (oil/water  Interface)  and  OA  (oil/air  interface).  A  positive  value  of  ST 
results  in  spreading. 

Gas  Chromatographic  Analysis 

To  estimate  the  fraction  of  the  oil  which  had  evaporated,  GC  analyses 
were  done  and  the  traces  compared  with  samples  which  had  been  evaporated  to 
known  extents. 
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A  0.5  wL  sample  was  Injected  into  a  Hewlett-Packard  Model -5750  gas 
chromatograph.  The  column  was  temperature  programmed  from  50°C  to  280°C 
at  15  degrees/min  with  a  one-minute  post  injection  interval.  The  column 
used  was  a  10  foot  long,  1/8  inch  O.D.  stainless  steel  column,  coated  with 
lOt  SE30  ultraphase  on  chromosorb  P,  A/W,  DMCS,  mesh  60/80.  A  flame  ioniza¬ 
tion  detector  was  maintained  at  320°C  and  the  injection  port  was  held  at 
300°C.  The  carrier  gas  and  the  hydrogen  flow  rate  were  at  about  30  mL/min. 

The  air  flow  rate  was  500  mL/min.  A  Shlmadzu  C-R1A  Chromatopac  integrator 
recorder  was  used. 

The  fuel  oil  samples  which  were  much  less  volatile  were  analyzed  in  the 
same  manner  using  a  HP  700  Laboratory  Chromatograph. 

Laboratory  Evaporation  (Air  Bubbling) 

To  prepare  oil  samples  of  known  extents  of  evaporation,  oils  were 
subjected  to  air  bubbling  under  controlled  flow  and  temperature  conditions. 
This  gives  data  which  can  be  used  to  estimate  environmental  weathering. 

In  addition,  samples  of  this  "artificially  weathered"  oil  were  taken 
and  selected  properties  determined;  namely,  the  distillation  curves  and  other 
property  curves  discussed  earlier. 

The  oils  were  air  bubbled  In  a  250  mL  graduated  cyclinder.  The  small 
cylinder  size  forced  the  use  of  low  air  flowrates,  generally  in  the  range  of  1 
L/min.  The  oil  was  evaporated  for  selected  "exposures"  to  obtain  evenly 
spaced  samples. 

To  take  each  sample,  the  cylinder  mass  and  the  volume  of  oil  at  both 
0°C  and  20°C  were  measured  and  approximately  40  mL  was  poured  into  the 
pour  point  apparatus.  This  40  mL  was  used  for  all  the  remaining  tests. 

The  cylinder  was  then  reweighed  and  the  volume  of  oil  was  measured  at  20°C 
only.  At  this  point  the  amount  of  air  which  passed  through  the  oil  was 
recorded. 

In  the  case  of  the  home  heating  oil  there  was  an  accidental  loss  of 
approximately  60  mL  during  the  experiment.  Fortunately,  there  was  still 
sufficient  oil  remaining  to  complete  the  experiment. 

Distillation 

The  distillation  method  used  (which  gives  oil  volatility  data)  is 
not  the  ASTM  method.  The  process  used  was  altered  in  order  to  predict  the 
evaporation  rate  more  accurately.  This  Is  done  by  measuring  the  liquid's 
boiling  temperature  not  the  lower  condensing  vapor  temperature.  The  flask 
from  which  the  oil  was  boiled  was  insulated  to  reduce  superheating,  and  the 
liquid  contained  glass  beads  (or  boiling  chips)  and  a  magnetic  stirrer  to 
maintain  uniformity.  A  thermocouple  was  used  to  measure  the  liquid 
temperature. 

The  amount  of  condensing  liquid  was  measured  against  boiling  tempera¬ 
ture.  Precautions  were  taken  to  prevent  the  evaporation  of  the  liquid  once 
It  reached  the  end  of  the  distillation  column. 


Flash  and  Fire  Points 


Flash  and  fire  points  were  measured  for  selected  oil  samples  to  provide 
an  indication  of  the  ease  of  ignition.  The  standard  ASTM  Cleveland  Open  Cup 
method  was  used  and  the  results  reported  in  °c. 

Emulsion  Stability 

The  wax  and  asphaltene  contents  of  the  fresh  oils  were  also  measured  as 
mass  fractions  by  standard  methods.  Assuming  that  these  components  do  not 
evaporate;  the  values  for  evaporated  oils  could  be  calculated  readily.  These 
quantities  are  used  later  in  estimating  stability  of  water-ln-oil  emulsions  as 
described  by  Mackay  and  Zagorski  (1982). 

In  addition,  the  fresh  oils  were  subjected  to  a  recently  developed 
emulsion  formation  and  stability  test  described  by  Mackay  and  Zagorski  (1982). 

In  total,  this  work  gave  a  comprehensive  set  of  data  which  could  be  used 
to  develop  expressions  relating  the  properties  as  discussed  in  Section  5.0. 


4.0  RESULTS 


4.1  Weathering  Experiments 

The  time  of  year  during  which  the  samples  were  weathered,  the 
condition  of  the  underlying  surface,  and  the  prevailing  weather  conditions 
during  certain  weathering  periods  have  led  to  specific  groupings  of  the 
samples  for  the  various  types  of  oil  as  outlined  below: 


Group 

Weathering  Period 

Sampl es 

: 

Prudhoe  Bay  Crude 

Winter  1980 

4  Feb  -  22  Feb  1980 

1  -  8 

Winter  1981 

1  Dec  -  4  Feb  1980 

9  -  18 

Winter  1981 
(emulsl fled) 

3  Feb  -  26  Feb  1981  and 

10  Mar  -  1.?  Apr  1981 

19  -  25  and 
34  -  35 

Spring  1981 

2  Mar  -  30  Apr  1981 

26  -  33 

No.  2  Home  Heating  Oil 

Winter  1980 

11  Feb  -  25  Feb  1900 

1  -  7 

Winter  1981 

6  Jan  -  25  Feb  1981 

8  -  18 

Spring  1981 

26  Feb  -  1 0  Mar  1 981 

19  -  30 

No.  4  Fuel  Oil 

Winter  1980 

18  Feb  -  6  Mar  1980 

1  -  6 

Winter  1981 

6  Jan  -  4  Feb  1981 

7  -  14 

Winter  1981 
(emulsified) 

3  Feb  -  26  Feb  1981  and 

2  Mar  -  6  Apr  1981 

15-21  and 
25,27,29,30 

Spring  1981 

2  Mar  -  13  Apr  1981 

22-24,26,28 

Tables  4.1  and  4.3  give  start  and  end  times,  duration,  average 
temperature,  wind  speed,  and  evaporative  exposure  (©)  for  each  sample,  e, 
the  evaporation  exposure.  Is  calculated  as  KAt/Vo  where 

K  Is  the  mass  transfer  coefficient  calculated  from  the  wind 
speed  (U(m/s))by  K  =  0.0025  U°*78  (Mackay  and  Matsugu,  1973) 

A  Is  the  oil  area  (m8) 

t  is  time  (s) 

Vo  is  the  Initial  volume  of  spilled  oil  (m8) 

The  first  (Winter  1980)  group  of  samples  were  weathered  during  cold 
weather  in  February-March  1980  when  the  ice  blocks  were  generally  Intact  and 
air  temperatures  were  generally  below  freezing.  The  samples  were  covered  with 
the  transparent  covers  to  protect  them  from  precipitation.  These  samples  were 
initially  deposited  in  cavities  in  the  ice  so  that  the  spill  area  was 
generally  less  than  the  box  area  for  the  first  few  days.  Although  these 
samples  started  out  on  a  solid  ice  surface,  by  the  end  of  the  period  some 
melting  had  occurred  so  that  a  simulated  melt  pond  situation  existed. 

The  second  (Winter  1981)  group  contains  those  samples  which  were 
weathered  during  very  cold  weather  in  January  1981  for  the  most  part  on  a 
solid  ice  surface.  They  were  not  initially  deposited  in  a  cavity  so  that  the 
oil  covered  the  major  portion  of  the  box  area.  The  samples  were  not  sheltered 
by  the  plastic  screens  so  that  snow  was  allowed  to  accumulate  on  the  samples. 
This  led  to  some  interesting  observations  of  oil/snow/ice  interaction  for  the 
darker  oils  (No.  4  and  Prudhoe  Bay  crude)  when  solar  radiation  and  warmer 
temperatures  melted  the  snow.  It  appears  that  the  darker  oils  are  initially 
sheltered  from  the  atmosphere  by  the  snow.  As  solar  radiation  warms  the  saw 
and  the  oil  underneath,  the  oil  appears  to  pocket  in  small  depressions  in  tne 
snow  and  ice  (up  to  a  few  cm  in  diameter)  (Figures  4.1  and  4.3).  It  is  quite 
possible  that  this  initial  sheltering  followed  by  the  collection  of  the  oil 
into  small  pockets  with  a  local  increase  in  slick  thickness  may  have  a 
significant  short-term  effect  on  the  weathering  rate.  The  lighter  No.  2 
home  heating  oil  did  not  show  this  pocketing  tendency.  It  saturated  the 
snow  evenly  shortly  after  being  covered. 

The  third  group  of  samples  (Winter  1981  emulsified),  refer 
primarily  to  those  samples  that  were  weathered  during  February  1981.  These 
samples  were  deposited  in  30  cm  x  30  cm  cavities  in  the  ice  to  prevent  runoff 
during  a  warming  period  in  February  1981.  The  samples  were  not  covered  with 
the  transparent  screens.  Melt  pond  conditions  prevailed  throughout  the 
period.  During  the  night  of  7-8  February,  heavy  rains  fell  inundating  the 
samples.  Although  the  samples  remained  intact  with  little  oil  loss,  some  oil 
emulsification  was  immediately  noticed.  The  emulsification  significantly 
affected  the  physical  properties  and  possibly  the  weathering  rates  of  the 
samples.  Accordingly,  some  of  the  later  samples  (March-Aprll  1981)  which 
showed  indications  of  emulsification  are  also  included  in  this  group. 


13 


Table  4.1  Exposure  Conditions  - 


Prudhoe  Bay  Crude. 


#  SAMPLE  START  TIME/DATE 


01 

PBC  *01 

1000 

04  FEB  80 

82 

PBC  *02 

100$ 

04  FEB  80 

03 

PBC  *03 

1000 

04  FEB  80 

04 

PBC  *04 

1000 

04  FEB  80 

OS 

PBC  *0S 

1000 

04  FEB  80 

06 

PBC  *06 

1000 

04  FEB  80 

0? 

PBC  *07 

1200 

20  FEB  80 

08 

PBC  *08 

1200 

20  FEB  80 

0? 

PBC  *09 

1100 

01  DEC  80 

10 

PBC  *10 

1100 

01  DEC  80 

11 

PBC  *11 

1400 

06  JAN  81 

12 

PBC  *12 

1400 

06  JAN  81 

13 

PBC  *13 

1400 

06  JAN  81 

14 

PBC  *14 

1400 

06  JAN  81 

15 

PBC  *15 

1400 

06  JAN  81 

16 

PBC  *16 

1400 

06  JAN  31 

17 

PBC  *17 

1400 

06  JAN  81 

18 

PBC  *18 

1400 

06  JAN  81 

19 

PBC  *19 

0900 

03  FEB  81 

20 

PBC  *20 

0900 

03  FEB  81 

21 

PBC  *21 

0900 

03  FEB  81 

22 

PBC  *22 

0900 

03  FEB  81 

23 

PBC  *23 

0900 

03  FEB  81 

24 

PBC  *24 

0900 

03  FEB  81 

25 

PBC  *25 

0900 

03  FEB  81 

26 

PBC  *26 

0900 

02  MAR  81 

27 

PBC  *27 

0900 

02  MAR  81 

28 

PBC  *28 

0900 

02  MAR  31 

29 

PBC  *29 

0900 

02  MAR  81 

30 

PBC  *30 

0900 

02  MAR  81 

31 

PBC  *31 

1100 

10  MAR  81 

32 

PBC  *32 

1100 

10  MAR  81 

33 

PBC  *33 

1100 

10  MAR  81 

34 

PBC  *34 

1100 

10  MAR  81 

35 

PBC  *35 

1100 

10  MAR  81 

36 

PBC  ST1 

— 

37 

PBC  ST2 

— 

38 

PBC  ST3 

— 

END  TIHE/DATE 

HOURS 

AVERAGE 

TEMP 

DEG  C 

AVERAGE  EVAP 

HIND  EXPOSURE 
M/SEC  X  iO"6 

0900  05  FEB  80 

024 

•1.2 

3.0 

.449 

0900  06  FEB  30 

048 

•5.2 

2.9 

.337 

0900  08  FEB  80 

096 

-3.0 

3.3 

1.947 

3.138 

0900  11  FEB  80 

170 

2.2 

3.0 

0900  19  FEB  80 

360 

-1.7 

3.4 

7.335 

0900  03  MAR  80 

672 

1.4 

3.2 

13.071 

0900  21  FEB  80 

021 

3.2 

i.l 

.185 

1400  22  FEB  80 

(50 

4.0 

2.1 

.707 

1000  02  DFC  80 

023 

-13.4 

2.1 

.320 

1000  03  DEC  80 

047 

-13  3 

4.0 

1.068 

1300  08  JAN  81 

047 

-3.3 

4.1 

1.134 

0800  10  JAN  81 

090 

•5.9 

3.0 

1.645 

0900  12  JAN  81 

139 

-8.1 

3.2 

2.730 

0900  15  JAN  Ai 

211 

-3.6 

2.8 

3.673 

1000  19  JAN  81 

308 

-7.5 

2.8 

5.450 

0900  22  JAN  81 

379 

-6.4 

2.6 

6.325 

0900  29  JAN  81 

547 

-4.3 

2.6 

9.090 

0900  04  FEB  81 

691 

-4.0 

2.8 

12.072 

0900  05  FEB  81 

048 

-B.2 

3.0 

.896 

0900  09  FEB  81 

144 

3.7 

2.9 

2.631 

0900  12  FEB  81 

216 

-1.6 

3.7 

4.685 

0900  17  FEB  31 

336 

1.6 

3.2 

6.494 

0900  19  FFB  81 

3B4 

-.7 

3.1 

7.189 

0900  23  FEB  81 

480 

1.0 

3.3 

9.118 

0900  26  FEB  81 

552 

1.6 

3.3 

10.926 

0900  04  MAR  Bi 

D4B 

-.4 

3.4 

i  Oil 

0900  09  MAR  81 

168 

1.0 

2.9 

3.042 

0900  12  MAR  31 

240 

1.5 

2.8 

4.230 

0900  16  MAR  81 

336 

1.8 

3.5 

7.071 

0900  19  MAR  81 

408 

1.2 

3.7 

3.329 

0900  23  MAR  81 

310 

1.0 

3.8 

6.895 

0700  26  MAR  31 

380 

1.4 

3.4 

7.631 

0700  30  MAR  81 

476 

2.3 

3.4 

9.713 

0800  06  APR  81 

645 

4.3 

3.5 

13.605 

0700  13  APR  81 

811 

5.1 

3.6 

17.443 

000 

— 

— 

0.000 

— “ 

000 

— 

— 

0.000 

non 


Table  4.2  Exposure  Conditions  -  No.  2  Home  Heating  Oil. 


♦  SAMPLE  START  TIME/OATE  END  TIME/DATE 


HOURS  AVERAGE  AVERAGE  EVAP 

TEMP  KIND  EXPOSURE 
DEG  C  M/SEC  X  iO"6 


01 

No2 

#01 

1000 

11 

FEB 

80 

0900 

12 

FEB 

30 

023 

-.2 

4.2 

.559 

0? 

No2 

*02 

1000 

11 

FEB 

eo 

1000 

13 

FFB 

80 

048 

-.8 

3.5 

1.021 

03 

No2 

#03 

1000 

11 

FEB 

80 

1000 

15 

FEB 

30 

101 

-3.2 

3.5 

2.105 

04 

No  2 

#04 

1000 

11 

FEB 

80 

1400 

18 

FEB 

80 

172 

-1.4 

3.5 

3.582 

05 

No2 

*05 

1000 

11 

FEB 

80 

0900 

25 

FEB 

30 

335 

.7 

3.0 

6.219 

06 

No2 

*06 

1000 

11 

FEB 

80 

SAMPLE  LOST 

000 

— 

— 

0.000 

07 

No2 

#07 

1200 

20 

FEB 

80 

1400 

21 

FEB 

80 

026 

4.1 

1.5 

.878 

oe 

No2 

#08 

1400 

06 

JAN 

81 

1300 

07 

JAN 

81 

,  023 

.4 

4.0 

.547 

09 

No2 

*09 

1400 

06 

JAN 

31 

0800 

10 

JAN 

81 

090 

-5.9 

3.0 

1.645 

10 

No2 

#10 

1400 

06 

JAN 

81 

0900 

12 

JAN 

81 

139 

-8.1 

3.2 

2.730 

11 

No2 

*11 

1400 

06 

JAN 

81 

0900 

15 

JAN 

81 

211 

-8.6 

2.8 

3.693 

1? 

No2 

*12 

140u 

06 

JAN 

81 

1000 

19 

JAN 

81 

308 

-7.5 

2.8 

5.450 

13 

No2 

*13 

1400 

06 

JAN 

31 

0900 

22 

JAN 

31 

379 

-6.4 

2.6 

6.325 

14 

No2 

*14 

1400 

06 

JAN 

81 

0900 

29 

JAN 

81 

547 

-4.3 

2.6 

9.090 

15 

No2 

#15 

1400 

06 

JAN 

81 

0900 

04 

FEB 

81 

691 

-4.0 

2.8 

12.072 

16 

No2 

*16 

0900 

03 

FEB 

81 

0900 

05 

FFB 

81 

048 

-8.2 

3.0 

.896 

17 

No2 

*17 

0900 

03 

FEB 

31 

0900 

09 

FEB 

81 

144 

-3.7 

2.9 

2.631 

J8 

No2 

:i8 

0800 

25 

FEB 

81 

1600 

25 

FEB 

81 

008 

5.0 

4.2 

.200 

19 

No2 

*19 

1100 

10 

MAR 

81 

0900 

11 

MAR 

81 

022 

2.9 

2.1 

.305 

20 

No2 

#20 

1100 

10 

MAR 

81 

0900 

12 

MAR 

81 

046 

2.4 

2.7 

.787 

21 

No2 

*21 

1100 

10 

MAR 

81 

1100 

14 

MAR 

81 

096 

3.2 

4.1 

2.850 

2? 

Ho2 

*22 

1100 

10 

MAR 

81 

1100 

16 

MAR 

81 

144 

2.6 

4.4 

3.633 

23 

No2 

#23 

1100 

10 

MAR 

81 

1200 

18 

MAR 

81 

193 

1.5 

4.5 

4.951 

24 

No2 

#24 

1100 

10 

MAR 

81 

0900 

19 

MAR 

81 

214 

1.1 

4.4 

5.366 

25 

No2 

#25 

1100 

10 

MAR 

31 

0700 

30 

MAR 

31 

276 

.8 

4.1 

6.566 

26 

No2 

*26 

1100 

10 

MAR 

81 

0700 

06 

APR 

81 

644 

4.3 

3.5 

13.584 

27 

No2 

*27 

1200 

16 

MAR 

81 

1100 

19 

MAR 

81 

071 

-2.1 

4.2 

1.739 

2B 

No2 

*28 

1200 

16 

MAR 

81 

0800 

23 

MAR 

81 

164 

-.5 

3.3 

3.235 

29 

No2 

*29 

1200 

16 

MAR 

31 

0800 

24 

MAR 

31 

188 

-.2 

3.0 

3.434 

30 

N#2 

*30 

1200 

16 

MAR 

81 

0800 

26 

MAR 

81 

236 

.6 

2.7 

4.016 

31 

No2 

ST1 

— 

— 

000 

— 

— 

0.000 

3? 

No2 

ST  2 

— 

— 

000 

— 

— 

0.000 

33 

No2 

ST3 

— 

— 

000 

— 

— 

0.000 
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Table  4.3 


Exposure  Conditions  -  No.  4  Fuel  Oil. 


» 

SAMPLE 

START 

TIHE/DATE 

END  TIHE/DATE  HOURS 

AVERAGE 

AVERAGE  EVAP 

TEHP 

HIND 

EXPOSUR 

DEG  C 

H/SEC 

X  10~4 

01 

No4  001 

1300 

18 

FEB 

80 

1300 

19  FEB 

80 

022 

.4 

5.3 

645 

02 

N#4  402 

1300 

18 

FEB 

80 

0900 

2(  FEB 

80 

042 

.8 

3.7 

.906 

03 

No4  403 

1300 

18 

FEB 

80 

1200 

22  FEB 

80 

(93 

2.8 

2.7 

1.579 

04 

No4  404 

1300 

18 

FEB 

80 

1400 

25  FEB 

80 

167 

3.0 

2.5 

2.677 

05 

No4  405 

1400 

03 

HAR 

80 

0900 

05  HAR 

30 

(43 

2.1 

3.0 

.784 

06 

No4  406 

1400 

03 

MAR 

80 

1500 

06  HAR 

80 

073 

3.2 

3.4 

1  493 

07 

No4  407 

1400 

06 

JAN 

81 

1300 

07  JAN 

81 

023 

.4 

4.0 

.847 

08 

No4  408 

1400 

06 

JAN 

81 

0800 

10  JAN 

81 

090 

-5.9 

3.0 

1.645 

0? 

No4  409 

1400 

06 

JAN 

81 

0900 

12  JAN 

81 

139 

•8.1 

3.2 

2.7.H 

10 

No4  410 

1400 

06 

JAN 

81 

0900 

15  JAN 

81 

211 

-8.6 

2.8 

3.693 

11 

Ho4  411 

14(0 

06 

JAN 

81 

1000 

19  JAN 

81 

308 

-7.5 

2.8 

5.4SP 

12 

No4  412 

1400 

06 

JAN 

81 

0900 

22  JAN 

81 

379 

-6.4 

2.6 

6.325 

13 

No4  413 

1400 

06 

JAN 

81 

0900 

29  JAN 

81 

547 

-4.3 

2.6 

9.090 

14 

H«4  414 

1400 

06 

JAN 

81 

0900 

04  FEB 

81 

691 

-4.0 

2.8 

12.072 

15 

No4  415 

0900 

03 

FEB 

81 

0900 

05  FEB 

81 

048 

-8.2 

3.0 

.396 

16 

No4  416 

0900 

03 

FEB 

81 

0900 

09  FEB 

81 

144 

-3.7 

2.9 

2.631 

17 

No4  417 

0900 

03 

FEB 

81 

0900 

12  FEB 

31 

216 

•1.6 

3.7 

4.685 

18 

Ni4  418 

0900 

03 

FEB 

81 

0900 

17  FEB 

81 

336 

-1.6 

3.2 

6.494 

1? 

No4  419 

0900 

03 

FEB 

81 

0900 

19  FEB 

81 

384 

-.7 

3.1 

7  139 

20 

Ho 4  420 

0900 

03 

FEB 

8i 

0900 

23  FEB 

81 

4B0 

1.0 

1.5 

v  418 

21 

No4  421 

0900 

03 

FEB 

81 

0900 

26  FEB 

81 

552 

1.6 

3.3 

10.926 

22 

Ho4  422 

0900 

02 

HAR 

81 

0800 

04  HAR 

81 

048 

-.4 

3.4 

1.011 

23 

No 4  423 

0900 

02 

HAR 

81 

0900 

09  MAR 

31 

16f 

1.0 

i  9 

3.042 

24 

Ho4  424 

0900 

02 

HAR 

81 

0900 

12  HAR 

81 

1.5 

2  8 

4.250 

25 

Ho4  425 

0900 

02 

HAR 

81 

0900 

16  HAR 

81 

336 

1.8 

3.5 

7.071 

26 

No4  426 

0900 

02 

HAR 

81 

0800 

19  HAR 

81 

408 

1.2 

3.7 

8.829 

27 

No4  427 

1100 

10 

HAR 

81 

0800 

23  HAR 

81 

310 

1.0 

3  8 

6.395 

28 

No4  428 

1100 

10 

HAR 

81 

0700 

26  PAR 

81 

380 

1.4 

3.4 

7.651 

29 

No4  429 

1100 

10 

HAR 

81 

0700 

30  HAR 

81 

476 

2.3 

3.4 

9.713 

30 

No4  430 

1100 

10 

HAR 

81 

0800 

06  APR 

81 

645 

4.3 

3  5 

13.605 

31 

No4  431 

1100 

10 

HAR 

81 

0700 

13  APR 

81 

811 

5.1 

3  6 

17.443 

32 

N«4  ST1 

— 

— 

000 

— 

— 

0.0(0 

33 

No4  ST2 

— 

— 

000 

— 

— 

0  000 

34 

No 4  ST3 

— 

— 

000 

— 

— 

0.000 

If 


Figure  4.3  Photograph  showing  closeup  of  dark  oil  sample  with 
oil  pockets.  Scale  is  in  centimeters. 
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The  fourth  group  (Spring  1981)  consists  of  those  samples  which  were 
weathered  during  March-Aprll  1981  when  wanner  weather  prevailed  (daytime 
temperatures  well  above  0°C),  and  clear  skies  led  to  significant  solar 
radiation  levels  at  the  test  site.  As  periods  of  melting  and  rainfall  were 
expected,  the  samples  were  deposited  in  the  30  cm  x  30  cm  cavities  and 
covered  with  the  transparent  screens.  These  samples  generally  experienced  a 
progression  of  underlying  surface  conditions  as  the  Ice  blocks  deteriorated. 
The  samples  would  rest  on  solid  ice  for  the  first  day  or  so,  spread  out  as 
melt  water  created  a  melt  pool  situation  after  a  few  days,  and  eventually 
weather  in  open  water  about  20-30  cm  deep  as  the  Ice  block  deteriorated 
completely.  This  hopefully  simulated  late  summer  conditions  on  the  North 
Slope  when  accelerated  weathering  of  a  spill  Is  likely  to  occur. 

The  No.  6  fuel  oil  samples  were  not  separated  into  groups  and  were 
analyzed  In  only  a  cursory  fashion,  as  very  little  weathering  was  observed. 

In  fact,  the  No.  6  oil  congealed  into  a  semi-solid  mass  at  the  colder  ambient 
temperatures  on  the  roof  and  did  not  spread  into  a  slick  as  did  the  other 
oils.  When  periods  of  warmer  temperatures  and  increased  solar  radiation 
melted  the  ice  around  the  oil,  the  semi -sol id  mass  would  become  partially 
submerged  in  the  melt  water.  If  the  melt  water  refroze,  the  oil  mass  would 
then  be  largely  encapsulated  in  the  ice,  further  sheltering  the  oil  from  the 
atmosphere.  The  fraction  of  oil  remaining  and  density  values  for  these 
samples  confirmed  the  absence  of  any  appreciable  weathering. 

4.2  Physical /Chemical  Properties 

The  results  of  the  physical  properties  analysis  conducted  for 
selected  oil  samples  by  the  University  of  Toronto  are  given  in  Tables  4.4, 

4.5,  and  4.6  for  the  Prudhoe  Bay  crude  oil,  the  home  heating  oil,  and  the 
heavy  fuel  oil.  The  distillation  curves  are  presented  in  Figure  4.4.  Gas 
chromatographic  traces  are  given  in  Figures  4.5  to  4.10  for  the  fresh  oil 
samples  and  for  selected  degrees  of  evaporation.  Figure  4.11  illustrates  the 
results  obtained  in  the  "purge  and  trap'1  analysis  for  water  solubility.  The 
emulsion  stability  test  results  are  given  in  Table  4.7. 

These  results  form  the  basis  of  the  model  calibration  work  described 
In  Section  5.0.  A  more  detailed  discussion  of  the  changes  in  physical  and 
chemical  properties  with  weathering,  and  the  implications  on  oil  spill 
countermeasures  and  cleanup  is  contained  in  the  report  by  Tebeau,  Meehan,  and 
Myers  (1982). 


2  Home  Heating  Oil — Fresh 


Figure  4.9 


Table  4.7.  Emulsion  Stability  Test  Results. 


OIL 

Test 

P.B.C. 

No.  2  H.H. 

No.  4  F.O. 

Mass  fraction 

a)  asphaltenes 

0.024 

0.0036 

0.032 

b)  waxes 

0.045 

0.029 

0.055 

12  hour  shaking 
test 

(1:3  oil-water 
ratio,  65  rpm 

oil  fraction 
emulsified 

0.5 

0.0 

0.8 

no  water 
layer 

(2  phases  only) 

forms  fairly 
stable  emulsion 

forms  no 
emulsion 

forms  stable 
emulsion 

5.0  MODEL  CALIBRATION  AND  DISCUSSION 


Having  obtained  the  analytical  and  property  data  for  the  oils  and  having 
developed  the  model  for  these  oils,  the  final  stage  was  to  parameterize  and 
run  the  model  equations  for  the  outdoor  test  conditions  at  Groton  and  predict 
the  rates  of  oil  evaporation  and  the  resultant  changes  In  oil  properties. 
Comparison  was  then  made  against  the  actual  properties  measured  experimentally. 
Because  the  outdoor  conditions  were  not  of  constant  wlndspeed  and  temperature, 
there  is  Inevitably  some  loss  of  fidelity  between  model  and  experiment,  and 
accommodation  is  possible  by  adjusting  the  rate  constant  for  evaporation  and 
the  environmental  conditions. 

In  practice,  the  full  oil  spill  model  was  not  used  because  spreading  and 
dispersion  did  not  occur.  Some  mousse  formation  did  occur  but  it  was  believed 
to  be  due  to  rainfall  and  is  not  properly  described  by  the  model  equations 
which  apply  to  oceanic  conditions.  Accordingly,  those  sections  of  the  model 
which  apply  to  weathering  and  physical  properties  were  removed  and  set  up  as  a 
separate  model . 

From  previous  work,  equations  were  proposed  for  the  dependence  of  the 
properties  on  temppr*ature  and  volume  fraction  evaporated  (F),  given  In  Table 
5.1.  These  equations  were  then  calibrated  to  the  specific  oil  types  by 
determining  the  equation  constants,  which  are  listed  In  Table  5.2.  These 
constants  were  derived  by  fitting  these  equations  to  the  experimental  data 
on  the  samples  obtained  from  the  air-bubbling  laboratory  weathering  procedures 
described  In  Section  4.1.  The  laboratory  weathering  data  were  used  to  derive 
the  constants  as  this  simulated  weathering  Is  caused  exclusively  by  evapora¬ 
tion,  which  Is  the  only  process  accounted  for  by  the  model  equations.  In  the 
Interests  of  simplicity,  an  attempt  was  made  to  fit  the  same  constants  (e.g., 

Kj ,  K?,  etc.)  to  all  three  oils,  but  this  was  not  possible  for  viscosity 
and  flash  and  fire  points. 

The  first  step  in  the  calibration  process  Is  to  relate  the  volume  frac¬ 
tion  evaporated  (F)  to  the  oil's  degree  of  evaporation,  1  .e. ,  the  evaporative 
exposure. 

The  evaporation  calculation  is  based  on  recent  work  by  Stiver  (1982)  and 
Mackay  and  Stiver  (1982)  and  is  an  extension  of  the  evaporative  exposure 
concept  developed  by  Mackay  and  Paterson  (1980). 

The  evaporation  rate  N  (m3/s)  is  given  by 
N  =  KAPv/RT 

where  K  Is  a  mass  transfer  coefficient  (m/s),  A  Is  the  oil  area  (m2),  P  Is 
the  oil  vapor  pressure  (atm),  v  Is  the  oil  liquid  molar  volume  (m3/mol),  R 
Is  the  gas  constant  (82  x  10"6  atm  m3/mol  K)  and  T  Is  absolute  temperature 
(K).  The  dimensionless  group  (Pv/RT)  Is  the  oil's  Henry's  law  constant,  l.e., 
the  ratio  of  the  concentration  In  air  (P/RT)  to  concentration  In  the  oil  (v), 
both  In  units  of  mol/m3.  This  group  varies  both  with  temperature  and  with 
extent  or  volume  fraction  of  evaporation,  F. 


Table  5*1  Properties  Measured  and  Correlated  against 
F  (volume  fraction  evaporated)  and  T  (tem¬ 
perature  K).  A  subscript  0  refers  to  conditions 
at  298  K  ie. 25  C  and  F  ■  0,ie. fresh  oil. 


PROPERTY 

UNITS 

FUNCTION 

Density 

Kg/m3 

DEN  -  DENq  (1  -  KX(T  -  T0) )  (1  +  K  F) 

Viscosity 

mPas 

(cp) 

VIS  -  VIS0  (exp  (K3(l/T  -  l/T0)exp(K4F) 

Aqueous 

Solubility 

g/m3 

SOL  -  SOLo  exP  (-K.F) 

Pour  Point 

°C 

PP  -  PP0  (1  +  K6F) 

Flash  Point 

°C 

FLP  -  FLF0  (1  +  K?F) 

Fire  Point 

°c 

FIP  -  FIP0  (1  +  KgF) 

Oil  Water  Inter¬ 
facial  Tension 

(3yne/cm) 

STW  -  STWo  (1  +  K  F) 

Oil-Air  Inter¬ 
facial  Tension 

mN/m 

(dyne /cm) 

STA  -  STA0  (1  +  K1QF) 

Wax  Content 

mass  fraction 

WAX  -  WAXq/  (1  -  F) 

Asphaltene  Con¬ 
tent 

- 1 

mass  fraction 

ASP  -  ASPq / ( 1  -  F) 
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Constanc 


PBC 


HH 


FO 


K1 

0.0008 

0.0008 

0.0008 

K2 

0.18 

0.18 

0.18 

K3 

9000 

3000 

7000 

K4 

10.5 

1.6 

5.0 

KS 

12.0 

12.0 

12.0 

K6 

0.35 

0.35 

0.35 

K7 

3.7 

0.35 

1.4 

K8 

3.7 

0.35 

1.4 

K9 

2.0 

2.0 

2.0 

K10 

1.0 

1.0 

1.0 

DENq 

895 

825 

905 

VISq 

35.0 

4.0 

23.0 

SOLq 

29.2 

3.0 

6.5 

PP0 

-2.0 

-27.0 

-3.0 

FLp0 

70.0 

100.0 

80.0 

FIPq 

80.0 

110.0 

90.0 

stw0 

27.0 

26.0 

30.0 

STAq 

30.0 

26.0 

32.0 

WAXq 

0.045 

0.029 

0.055 

ASPo 

. 

0.024 

0.0036 

0.032 

If  the  oil's  Initial  volume  Is  V0(m3)  then 
N  *  V0dF/dt 

hence  dF/dt  -  KAPv/V0RT 

or  dF/dO  *  Pv/RT  where  0  ■  KAt/V0 

0  Is  a  dimensionless  group  termed  the  "evaporative  exposure." 

The  problem  Is  now  to  obtain  an  equation  relating  Pv/RT  to  F  and 
temperature.  The  simplest  procedure,  which  Is  suitable  for  fairly  volatile 
oils.  Is  to  obtain  the  atmospheric  pressure  distillation  equation  of  boiling 
point  vs.  F,  namely, 

TB  *  A,,  +  BnF  (Af,  and  Bn  are  coefficients) 

For  several  oils,  the  following  relationship  applies  to  give  (Pv/RT)  at 

22°C: 

In  (Pv/RT)  *  6.3  -  0.035  TB 

hence  In  (Pv/RT)  -  6.3  -  0.035A,,  -  0.035  BnF 

and  substituting  from  above 

In  (dF/dO  )  *  6.3  -  0.035  A„  -  0.035  BnF 

Integrating  this,  and  applying  the  Initial  condition  that  0  «  0  at  F  -  0,  gives 

0  *  (exp(0.035  B„F)  -  1 )/(0.035  Bn  exp  (6.3  -  0.035  A„)) 

which  gives  an  analytical  relationship  between  0  and  F,  the  evaporative 
exposure  and  mass  fraction  evaporated. 

This  approach  was  satisfactory  for  the  Prudhoe  Bay  crude  oil, but  It 
failed  for  the  distillates  because  of  their  high  boiling  points.  Accordingly, 
separate  constants  were  fitted  for  these  oils,  the  0.035  and  6.3  being  changed. 
The  final  constants  are  given  below. 

Prudhoe  Bay  crude  oil 

In (Pv/RT)  -  6.3  -  0.035TB 

Tb  »  A]  +  B]F  Aj  ■  439  Bi  ■  804 

Home  Heating  Oil 

In (Pv/RT)  ■  50.1  -  0.1 2TB 

Tfi  *  A2  +  B2F  A2  *  521  B2  *  126 

No.  4  Fuel  Oil 

In (Pv/RT)  ■  50.1  -  0.1 2TB 

TB  ■  A3  ♦  B3F  A3  »  518  B3  -  239 
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It  is  hoped  In  future  work  to  fit  a  single  equation  to  cover  the  entire 
range  of  oils  Tlkely  to  be  encountered  environmentally. 

For  temperatures  TE0O  other  than  22°C  (295°K),  a  correction  Is 
applied,  namely 

ln(PE/P22)  *  10.6  TB(l/295  -  1/TE) 

Here  TE  Is  the  actual  environmental  temperature. 

Now  v  also  changes  with  temperature  (as  density  p  changes)  as  follows: 

pe  =  P22  H  -  Kl<TE  -  295)) 

therefore  vE  =  v22/(1  -  Ki(Te  -  295)) 

Combining  these  two  corrections,  and  assuming  that  In  (RT)E  %  1  n (RD22 
yields 

ln(Pv/RT)E  =  ln(Pv/RT)22  -  ln(l  -  K] (TE  -  295))  -  10.6TB(l/295  -  1/TE) 


The  constant  K]  is  given  In  Table  5.2. 

In  this  manner,  a  set  of  calibrated  equations  was  obtained  for  each 
oil  type,  which  calculates  the  volume  fraction  evaporated  (F)  from  the 
"evaporative  exposure,"  and  then  calculates  the  oil  physical  properties 
as  a  function  of  F.  The  next  step  in  the  modelling  effort  was  to  compare 
the  predicted  values  of  F  and  the  oil  physical  properties  to  the  experimental 
data,  to  judge  the  accuracy  of  the  model  equations.  A  complete  tabulation  was 
prepared  (which  Is  not  presented  here,  but  is  available  at  the  RAD  Center)  of 
the  individual  experimental  and  calculated  values  for  each  oil  sample.  The 
summary  results  are  presented  in  Figures  5.1  and  5.19  in  the  form  of  graphs 
of  predicted  versus  experimental  property  values  for  volume  fraction 
remaining,  density,  viscosity,  solubility,  pour  point,  flash  point,  fire 
point,  and  interfacial  tensions.  The  data  for  each  oil  type  are  plotted 
separately. 

The  plots  for  volume  fraction  evaporated  are  given  in  Figures  5.1  to 
5.3.  The  predicted  values  were  obtained  by  first  computing  the  evaporative 
exposure  (e).from  the  time  duration,  average  wind  speed,  and  average  temper¬ 
ature  values  in  Tables  4.1  to  4.3,  and  then  using  evaporative  exposure  to 
compute  F  according  to  the  formulation  outlined  above.  Experimental  data,  for 
the  air-bubbled  samples.were  measured  directly  in  the  laboratory.  However, 
such  data  were  unavailable  for  the  samples  weathered  under  environmental 
conditions  at  the  RAD  Center  due  to  inaccuracies  In  the  measurement 
procedures.  For  these  samples,  volume  fraction  evaporated  could  only  be 
estimated  by  Inspecting  the  GC  traces  and  matching  with  samples  of  known  F. 

The  estimated  F  values  are  given  In  Tables  4.4  to  4.6. 

The  comparisons  between  predicted  and  calculated  oil  physical  properties 
are  given  In  Figures  5.4  to  5.19.  Predicted  values  were  calculated  according 
to  the  formulations  set  forth  In  Tables  5.1  and  5.2.  Experimental  values  were 
obtained  from  the  University  of  Toronto  analysis,  the  results  being  summarized 
In  Tables  4.4  to  4.6.  For  the  density  end  viscosity  plots,  values  for  samples 
which  were  subject  to  emulsification  have  been  Identified,  as  these  experi¬ 
mental  values  are  not  expected  to  agree  with  the  values  predicted  by  the 
evaporation  model. 


Experimental  Fv  Based  on  Exposure 


FIGURE  5.1  FRACTION  EVAPORATED  BY  VOLUME— PRUDHOE 


FIGURE  5.2  FRACTION  EVAPORATED  BY  VOLUME-NO.  2  HOME  HEATING  OIL 


I 
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FIGURE  5.3  FRACTION  EVAPORATED  BY  VOLUME-NO.  4  FUEL  OIL 


Predicted 

(g/cm3) 


FIGURE  5,5  DENSITY— NO. 2  HOME  HEATING  OIL 


Density  Measured  (g/cm^) 


FIGURE  5.6  DENSITY-NO.  4  FUEL  OIL 


Viscosity  Measured  (cp) 


FIGURE  5,7  VISCOSITY— PRUDHOE  BAY  CRUDE 


Viscosity 

Predicted 

(cp) 


FIGURE  5.8 


jsity  Measured  (cp) 


—NO.  2  HOME  HEATING  OIL 


10 


Viscosity 
Predicted 
(cp)  40 
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A  Indicates  Emulsified  Sample 


60  80  100 


Viscosity  Measured  (cp) 


FIGURE  5,9  VISCOSITY-NO.  FUEL  OIL 


FIGURE  5.10  AQUEOUS  SOLUB I L I TY — PRUDHOE  BAY  CRUDE 


Solubility 

Predicted 


□  AB 1530 


FIGURE  5.11  AQUEOUS  SOLUBILITY-NO.  2  HOME  HEATING  OIL 
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A87260,27 

12  V 

,  AB1770 
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*  /OAB667 
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18  24 
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Pour  Point  Measured  (3C) 
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FIGURE  5.13  POUR  POINT— PRUDHOE  BAY  CRUDE 


Pour  Point  Measured  (~C) 


FIGURE  5.15  POUR  POINT-NO.  4  FUEL  OIL 
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FIGURE  5,16  FLASH  AND  FIRE  POINT 


Interfacial  Tension  Measured  (N/m) 


FIGURE  5.17  INTERFACIAL  TENSION— PRUDHOE  BAY  CRUDE 


FIGURE  5,18  INTERFACIAL  TENSION-NO.  2  HOME  HEATING  OIL 
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Tension 

Predicted 
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FIGURE  5.19  INTERFACIAL  TENSION— NO.  4  FUEL  OIL 


The  agreement  between  the  predicted  values  and  those  measured  as 
Indicated  In  these  figures  was  judged  to  be  satisfactory  In  most  cases. 
Prediction  of  the  extent  of  evaporation  Is  reasonably  accurate  considering  the 
variability  in  environmental  conditions  and  the  uncertainty  of  the  measured 
extent  of  evaporation.  Density  predictions  were  satisfactory  with  the  excep¬ 
tion  of  a  few  samples  In  which  the  discrepancy  is  attributed  to  formation  of 
water-ln-oll -emulsions  which  formed  because  the  samples  were  exposed  to  rain. 
This  phenomena  Is  believed  to  cause  the  only  major  discrepancies  in  the  vis¬ 
cosity  predictions  as  well.  Aqueous  solubility,  pour  point,  and  flash  and 
fire  point  estimation  gave  good  results.  Interfacial  tensions  were  unpre¬ 
dictable.  This  is  probably  due  to  oxidation  occurring  In  the  environment 
which  was  not  reproduced  in  the  lab.  Under  the  prolonged  influence  of  sun¬ 
light  and  oxygen.  It  is  likely  that  certain  components  of  the  oil  (such  as  the 

asphaltenes)  oxidize  and  form  surface  active  compounds  which  reduce  the  Inter¬ 

facial  tensions. 

In  general,  the  model  provides  a  satisfactory  predictive  capability.  The 
accuracy  could  be  improved  by  selection  of  better  equations  and  incorporation 
of  more  experimental  data;  but  the  significant  advance  in  this  work  has  been 

to  establish  for  the  first  time  the  predictive  framework  and  obtain  reasonable 

parameter  values. 

6.0  CONCLUSIONS 

An  oil  spill  model  has  been  developed,  described,  and  provided  to  the 
Coast  Guard  RAD  Center,  Groton,  CT.  The  model  as  originally  developed 
Includes  calculation  of  rate  of  spreading,  drift,  evaporation,  dispersion  and 
water-ln-oll  emulsion  formation.  It  also  calculates  the  oil  physical  proper¬ 
ties  (density,  viscosity,  pour  point,  aqueous  solubility,  flash  point,  fire 
point,  and  Interfacial  tension)  as  a  function  of  these  weathering  processes. 

The  focus  of  this  study  was  to  develop  and  calibrate  a  simplified  version 
of  this  model  which  predicts  volume  fraction  evaporated  and  physical  properties 
based  on  evaporation  alone  for  three  types  of  oil  under  arctic  conditions.  The 
simplified  evaporation  model  was  calibrated  using  data  obtained  In  simulated 
weathering  (air-bubbling  evaporation)  experiments  in  the  laboratory.  The 
accuracy  of  the  model  was  then  judged  by  comparing  the  predicted  values  to 
experimental  data  obtained  from  outdoor  weathering  experiments  at  the  RAD 
Center  during  the  winters  of  1979/80  and  1980/81.  The  results  indicate  that 
the  model  predicts  oil  density,  viscosity,  pour  point,  and  aqueous  solubility 
with  a  reasonable  degree  of  accuracy.  Attempts  to  predict  flash  and  fire 
points,  and  Interfacial  tensions  were  not  satisfactory.  In  view  of  this, 
the  simplified  model  represents  a  significant  contribution  to  predicting  oil 
behavior  and  properties  under  arctic  marine  environmental  conditions. 


56 


REFERENCES 


Huang,  J.C.  and  Monastero,  F.C.,  1980.  Review  of  the  State  of  the  Art  of 
Oil  Spill  Simulation  Models,  unpublished  report  to  American  Petroleum 
Institute. 

Jordan,  R.E.  and  Payne,  J.R.,  1980.  "Fate  and  Weathering  of  Petroleum 
Spills  in  the  Marine  Environment,"  Ann  Arbor  Science  Publishers, 

Inc.,  Ann  Arbor. 

Kirstein,  B.  and  Payne,  J.R.,  1982.  Private  Communication,  SAI,  Ltd. 

Mackay,  0.  and  S.  Paterson,  1978.  "Oil  Spill  Modelling,"  Proceedings  of  a 
Workshop  held  in  Toronto,  Canada,  Pub.  No.  EE-12,  Institute  for 
Environmental  Studies,  University  of  Toronto. 

Mackay,  D.,  Buist,  I.,  Mascarenhas,  R. ,  and  Paterson,  S.,  1979.  "Oil 
Spill  Processes  and  Models,"  Environment  Canada  Environmental  Protection 
Service,  Publication  EE -8. 

Mackay,  D.  and  Hossain,  K.,  1982.  "Interfacial  Tensions  of  Oil,  Water, 
Chemical  Dispersant  Systems,"  C.J.  Ch.E. 

Mackay,  D.  and  Matsugu,  R. ,  1973.  "Evaporation  Rates  of  Liquid 
Hydrocarbon  Spills  on  Land  and  Water,"  J.  Chem.  Eng.,  Vol .  51. 

Mackay,  D.  and  Paterson,  S. ,  1980.  "A  Mathematical  Model  of  Oil  Spill 
Behaviour,"  Report  prepared  for  Environmental  Emergencies  Branch, 
Environmental  Protection  Service,  Publication  EE-7. 

Mackay,  D.,  Paterson,  S.  and  Nadeau,  S.,  1980.  "Calculation  of  the 
Evaporative  Rate  of  Volatile  Liquids,  "Proceeding  from  the  1980  National 
Conference  on  Control  of  Hazardous  Material  Spills. 

Mackay,  D.  and  Stiver,  W.,  "Evaporation  of  Crude  Oils"  (in  preparation). 

Mackay,  D.  and  Stiver,  W.,  1982.  "Artificial  Weathering  of  Crude  Oils," 
Spill  Technology  Newsletter,  Vol.  7. 

Mackay,  D.  and  Zagorski,  W.,  1982.  "Studies  of  Water-in-Oil  Emulsions," 
EE-34,  Environmental  Emergency  Branch  of  Environment  Canada. 

Stiver,  W.,  1982.  "Studies  of  Crude  Oil  Evaporation,"  B.A.Sc.  Thesis, 
University  of  Toronto. 

Tebeau,  P.  A.,  Meehan,  T.M.,  and  tyers,  J.  C.,  1982.  "A  Laboratory  Study 
on  Oil  Weathering  Under  Arctic  Conditions,"  U.  S.  Coast  Guard  R4D  Center, 
Groton,  CT. 


57 


APPENDIX  A 


OIL  SPILL  BEHAVIOR  MODEL 


INTRODUCTION 

This  Appendix  gives  a  brief  description  of  the  model  and 
Instructions  for  input  data.  The  program  could  be  modified  to  become 
interactive  by  insertion  of  suitable  statements. 

The  oil  parameters  used  here  are  purely  illustrative  and  do  not 
reflect  the  properties  of  a  given  oil. 

The  input  data  are  grouped  into: 

(a)  Computing  parameters 

(b)  Spill  quantities 

(c)  Environmental  conditions 

(d)  Oil  properties 

(e)  Process  constants 

(a)  Computing  Parameters 

The  following  quantities  must  be  defined: 

DTIM  -  time  increment  (usually  100  s) 

TCOMP  -  total  time  to  be  computed  (e.g.,  24  hours) 

POFH  -  print  out  frequency  in  hours  (e.g.,  every  0.1  hours) 

(b)  Spill  Quantities 

VOLTS  -  total  volume  spilled  In  nr*  (eg.  100) 

TIMSP  -  duration  of  spill  in  seconds.  This  should  be  a 
multiple  of  DTIM.  For  an  instantaneous  spill  set 
TIMSP  to  DTIM. 

TTK  -  initial  oil  slick  thickness  (m)  usually  set  at 
0.02  m  (2.0  cm) 

TTN  ■  -  thickness  of  sheen  (m)  usually  0.000005  m  (5  x  10"®  m) 

(c)  Environmental  Conditions 

TC  -  temperature  °C  (eg.  1 S°C ) .  Note  that  CF  Is 
printed  out  also. 

WSKH  -  wlndspeed  in  km/h  (eg.  20  km/h).  Note  that  speeds 
In  m/s  and  knots  are  also  printed  out. 

FDR  -  oil  drift  factor  (usually  0.035  or  3.5%  of  wind  speed) 


A-l 


Note  that  all  values  must  be  inserted  at  25<>C.  A  correction  Is 
then  made  to  the  ambient  temperature  for  density,  viscosity  and  vapor  pressure. 

FMAX  -  maximum  fraction  evaporated  (usually  0.4  but  possibly 
lower) 

ST  -  oil -water  Interfacial  tension  (usually  24  4ynes/cm) 

DEN25  -  fresh  oil  density  at  250c  (approximately  850  kg/m3) 

VIS25  -  fresh  oil  viscosity  at  25oc  (approximately  10  cPolse) 

SOLO  -  fresh  oil  solubility  (approximately  30  g/m3) 

SENT,  VIST  -  constants  relating  density  and  viscosity  to 
temperature 

DENK,  VISK,  SOLK  -  constants  relating  the  change  in  density, 

viscosity  and  solubility  to  weathering.  Typical  values 
are  100,  4.0  and  12  respectively.  The  values  of  density 
viscosity  and  solubility  (DENX,  VISX,  SOLX)  at  10% 
weathering  are  calculated  and  printed  out  as  a  check. 


(e)  Process  Constants 


Evaporation 


None  are  inputed  but  CE1 ,  the  evaporation  mass  transfer  coefficient 
Is  calculated  from  wind  speed.  It  Is  generally  about  0.005  to  0.01  m/s.  CE2 
and  CE3  are  A  and  B  respectively  taken  from  the  distillation  curves  for  each 
oil  (See  Figure  4.4). 

Dispersion 

Three  constants  are  Inserted,  CD1 ,  CD2,  and  CD3  with  values  0.00005, 
2.0  and  1.0  respectively. 

Emulsion  Formation 

Three  constants  are  Inserted,  CM1  with  value  of  0.65,  CM2  with  value 
of  2  X  10~6,  and  CM3  with  value  of  0.6  for  PBC,  0.25  for  HH,  and  0.7  for  FO. 

Spreading 

Four  constants  are  Inserted,  CS1 ,  CS2,  CS3,  and  CS4  with  values  1.0, 
150,  5.0,  and  8.0  respectively. 


Note 

In  some  cases,  the  values  of  these  constants  are  oil  dependent. 


Shut-Off  Factors 


Five  factors  are  read  In  with  a  value  of  1 .0  If  the  processes  of 
evaporation,  dlsperlson,  emulsion  formation,  spreading,  and  drifting  are  to 
proceed  as  normal.  A  value  of  zero  stops  the  process  and  thus  permits  the 
Impact  of  that  process  on  overall  spill  behavior  to  be  assessed. 


COMPUTING  PROCEDURE 

1  -  the  program  Initializes  a  number  of  variables  then  proceeds  to 

calculate  the  Initial  disposition  of  the  slick  and  the  oil 
properties. 

2  -  a  print-out  Is  then  made  of  the  spill  conditions 

3  -  the  oil  spill  behavior  Is  then  calculated  as  follows  for  the 

thick  slick  and  the  sheen 

4  -  time  change  Is  calculated 

5  -  evaporation  amounts  and  rates  are  calculated 

6  -  dispersion  amounts  and  rates  are  calculated 

7  -  emulsion  formation  amounts  and  viscosity  are  calculated 

8  -  spreading  Is  calculated 

9  -  the  new  oil  volumes,  areas,  composition  and  properties  are 

calculated,  any  "new"  spilled  oil  being  added 

10  -  the  slick  position  and  shape  Is  calculated 

11  -  the  condition  of  the  slick  Is  then  printed  out  If  requested 

12  -  the  program  then  returns  to  step  4,  checking  to 

determine  If  the  program  should  stop. 
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aARREN  STIVER  * 


CIL  SPILL  MCUEL 

NC.  2  HCME  ME AT INw  OIL 


SET  PAR AME  TCP  VALUES 

SET  COMPUTING  CCNCITICNS 

TIME  INCREMENT  IN  SECONDS  -  071* 

CT  I  M*  l  00  • 

TOTAL  TIMS  COMPUTED  IN  HCUP5  r  TCCMF 

TCCMP  »  23. C 

M7T«TCCMO*260C/DTI* 

PRINT  OUT  FPE3UENCY  IN  HOURS 

POPH  *  4.0 

NLP*PONH»2eCO/CT  IN 

MT*0.0 

LOCP-NLP 

SET  SPILL  CONDITIONS 

1BTAL  VOLUME  SPILLEC  -  VCLTS  (M3) 

VOLTS-  100 

Curat i on  cr  spill  in  seccncs  (min  is  one  increment)  -  timsp 

ti MSP— i oo . c 
TIMhP-T  IMSp/2tCC. 

TIM-0.0 

CVOL-VCLTS/TIMSP 
Vl*0VCL*0T  ltd 

INITIAL  THICK  SLICK  THICKNESS  -  7  T  K  (Ml 
TTK-0.C2 

SHEEN  THICKNESS  -  TTN  (M) 

TTN*0. 0O0SC5- 
20L-O.C 
200-0.0 
VRAT-0.0 

SET  ENVIRONMENTAL  CCNCITICNS 
TEMPERATURE  -  TC  (DEGREES  C> 

TC  »  2, 

TK*TC-*273. 

TF*32.0h1.€*TC 
V I NCSPEEO  -  V3KM  (KM/HR) 

VSKH-20 
4SMS-VSKH/  ;  .C 
»SKN**SMS •  1  .042 
OR  IPT  FAC  TCP  -  FOR 
FOR-O ,075 
vor-foc-asm* 

THETA  «  0. 


2  t 

2  S 

2  C 
31 
22 
22 

34 

2  ? 
2t 

37 

3  fc 

35 

4  C 

4  1 
42 
42 
44 


4 1 
47 
4( 
4S 


SC 
S  1 


s 

9t 


e  i 


0  2 
22 
C  4 

I! 

(( 


<7 
66 
6< 
7  C 


c 

c 

c 

c 


c 


c 


c 


c 


c 

c 


c 

c 


c 

c 


c 

c 


c 


SET  CIL  PROPERTIES 

MOLAR  VOLUME  -  VMdL  tM3/*CL> 

MAX  FRACTION  2VAP3RATSC  -  F-AX 

ahp  *  o.i2  -  ic.e  •  <  i.c/2<ss.  -  i.c/tki 

EHP  -  50.1  -  AUOGC  1.0  -  C.00C24ITK  -  5SS.JI 

FMAX  ■  0  •  J 

FTN-O.C 

FT K* 0.0 

Pa  0. 000002 

CIL  MATER  INTEPFACIAL  TENSION  -  «T  CCVNES/CM) 

STa  26.0 
STK-ST 

stn*st 

CIL  DENSITY  -  CENC  ( KG/M  3  ) .  aEATnERING  CONSTANT  -  DENK 
CEN2S  •  62  C • 

cent-  o.occe 

CeNO  ■  0EN23  •  <1  -  CENT 4<TC-25  )  ) 

CE  NK  ■  0  •  l  6 

CIL  VISCOSITY  -  Vise  (C?0ise>.  MIATFERING  CONSTANT  -  VISA 
V1S2S  *  3.4 
VIST  *  3  0-JC. 

VI  SK  •  1  .6 

Vise-  VISJ*  *EXP( vl ST*(  1 ,0/TK  -l. C/2S3.  I) 

CIL  SULUfl  lL  I T Y  -  jOLl  (G/M2I.  *£ATh£RING  CONSTANT  -  SOLK 
SOLO  -  3.0 
SOLK  «  1 2 • C 

CE  NX  -  OENC  •  <1.  ♦  CENK  A  0  •  1  > 

VISA-  VISO  -  EXP  (VISK-O.t) 

SOLX  A  SOLC  •  EXPI-  SOLK*  3.1) 


8V  APORAT I C  N  CONSTANTS 
CEl-0.000e7*(ASKh)AAC.76 
CC2-  321 • 

CE 3  a  126. 

VTKE-0.0 

VTNEaO.O 

VTE-O.C 


Cl SPERSI CN 
COl-O.OCJ  )! 
C02-2.0 
CD  3* 1 .0 
VTKCaO .0 
vTNOaO.O 
VTQaO.O 


CONSTANTS 


EMULSICN  CCNSTAAT3 

CMlAO.eS 

0*2  a  2. 02-06 

CM3  -  0 . 2  S 

Am  aO  .0 

OF  MAQ.C 


SPREAOING 
CSl«l .0 
CS  2*130. 
CS3*3.0 
CS4*8.0 


CONSTANT, 


A-12 


VOLTS  IS  TCTAL  VOLUME  SPlLLcC  (CUBIC  KCTREtl 

VS PL  IS  VCLLMu  SP ILL  cO  TC  GIVEN  IlME  (CUBIC  METRES) 

Tt  M  *  TIMM.  ItMM.TINC  ACE  T  *  ME  (SCCCNCS*  RINLTES*  NOUNS*  OATS  I 
CTIM  IS  TIME  INCREMENT  (  SECCNOl  I 

Tt  PSP  IS  -SPILL  DLnATIQS  (SSCCkCS)  (NLST  iucceo  txmi 
CVOL  13  SPILL  VCLCPC  INCACMENT 

PMAX  IS  MAXIMUM  FRACTION  CP  Tm€  CIL  *M  ICH  CAM  EVAPORATE 
•SKH  •  »SMS  ANO  «SiCN  Add  m  IN-  SPEEDS  IN  MM/P.  M/S  AMO  MNOTS 
FOR  IS  MlNC  cmipt  pactlr 

TTn.  TTM.  T-C  arc  ASE  thin,  ">IC«  anc  AVERAGE  spill  thicknesses 
ATM.  ATM.  <  TO  ARC  THIN.  TflCN,  ARC  TCTAL  SP  ILL  AREAS 
ISOLARE  METRES) 

VTN.VTK,  VTC  ARfc  TM IN*  TUCK  ARC  TCTAL  SPILL  VOLUMES 
(CUBIC  METRES) 

PTn,  FT*  AAC  ft  ARfe  THIN,  THICK  ANC  TCTAL  FRACT ICNS  tVAPCPATEO 
ST,  STN.  AAC  STM  ARC  INITIAL*  *FC  TP  IN  ANC  THICK  SL  i  CK  INT2A FACIAL 

oInc!°VISC  AND  SCLO  ARE  INITIAL  CENSITV,  VISCOSITY.  AND  SOLUBILITY 

CENT?CV?ST.  VPRT  ARE  CCNSTANTS  RELATING  DENSITY.  VISCOSITY  ANC 
VAPOH  PRESSLRb  TC  TEMPERATURE 

CCN2S,  VtSif*  VPH2S  ARE  FRESH  CIL  DENSITY*  VISCOSITY  AND 
VAPCP  PRESSURE  AT  25  DEGREES  C 

Cfc NX*  VISA.  SOLA,  VPRX  ARC  VALLES  CP  DENSITY.  VISCOSITY  AND 

VARCR  PRESJLPG  AT  IQS  VEAtHER |NG 

INITIAL  RAT  INDICATES  PATES  CP  PROCESSES 

2L  •  ZW  ARE  LENGTH  AND  NlOTM  CP  SLICK 

ZOL  IS  DRIFT  LENGTH  CFTHE  SPILL  SCLRCE  PClNT 

2DO  IS  OR  I  AT  LENGTH  FRCP  LAST  SPILL  PCINT 

TEXP  IS  EXPOSURE  TIME  ( SECCNOS ) 

INITIAL  D  MEANS  CELT  A  CR  CIPFEREPCE 
PINAL  E  INDICATES  EVAPORATION 
PINAL  C  INDICATES  OtSPERSICN 
PINAL  M  INDICATES  EAcLSICN 
PINAL  S  INDICATES  SPREADING 
PINAL  V  INDICATES  OR  IF  T 

SHUT  OFF  FACTORS 
SOP PC*  1.0 
SOPFDm  l  .0 
SOFFM»i .0 
50FF5-1 .0 
SOP F*»*  I  ,0 

SHUT  OFF  PROCESSES 

PMAX»FMAX*SCFFE 

CEi*SOFFe»CE l 

COl»SOFFO*CCl 

CM2»SCFFM»CH2 

C$  t«CS !■ SCP  A  S 

CS2-CS2*S0AFS 

por»forpscpf» 
initial  CQNClTtCNS 

CALCULATION  OF  INITIAL  DISPOSIUCN  CP  01LSL1CK 
VSPL»VI 

ATK-VI/(TTPVC54*TTN) 


ATN>ATK*C34 

ATQmATKAATN 

VT*-ATk*TT« 

VTN*ATN*TTN 

C 

C  CALCULATE  INITIAL  OIL  PROPERTIES 

C  CENSITV 

CCN  *  CENO«(I.  ♦  CENP-FTM) 

C  VISCOSITY 

VIS  •  VISC  •  EXP ( V I SK*PTK  ) 

VI  SEN* VIS 
C  SOLLOIL  ITV 

SOL»SOLC-d  XP ( -SCLKMPTK  I 

VRlTE(e*SOS)  VCLTS  L 

SOS  FORMAT  I • 1 •• 'TOTAL  VOLUME  TC  BE  IP ILLCO* • .Ft .2 *2X.**CU0IC  METRES*/// 
l  /✓) 

NRITEI6.7I2)  TIMSP.TIMHP 

712  format ( •  •••duration  of  spill  (s  )  • • ijx.fo.i.sx. *or *.pi i • j.ix.'h* > 
NR  1TE( A  •  7  t  2  )  VK 

7|3  FORMAT (•  •••SPILL  VOLUME  INCREMENTS  CMS  )  •  •  5 A *P«.  I ) 

*R!TE(6*M7>  T7M*TTN,C'i% 

7»7  FORMAT  (  •  •••PIITHL  5PILL  THICKNESSES  (  M )  »  ,SX  ,P  |  i  .«  ,2X  ,  •  TH  ICR 
•SLICK* *SX*Fll ,C*2A* 'THIN  SL I CK •  •  3* • P  1 1 . 6 • 2X . • ARC A  FACTOR*//) 

VR I  TE  (  ft  •  7  |  •  ) 

715  FORMAT ( *0*  •  •ENVIRONMENTAL  CONO 1 1 ICN f  •/) 

BR|TE(6,35C>  TC.Tr 

35  C  FORMAT (• C •••  TdMPEAATLRE  (DEGREES  C)*«t0X.P9*t ,SX  **QR* .PB  .1 .2 A . 

•  •(DEGREES  PI*  ) 

VP[TE(6  *3S I  )MSKH  ,«SMS*MSKN 

JSl  FORMAT ( •  •••VINO  SPEED  ( KM/H ) *  .  1 1 A  ,P * . I  • 5X  •  • CR* .P«. I .2* • • M/S • . SX . 
•*OR*.P«.I . I X  * • KNOTS •  I 
VRITE(«*7IC)  PCR 

710  FORMAT  (  •  •.'DRIFT  FACTOR*  *2CX«PE.3//I 
VRITE(b.7QC) 

7  CO  FORMAT ( *0*  #  '  C I L  PROPERTIES*/) 

NRITE(6.P01) 

701  FORMAT! *0*  *33X  •• INI T  1AL  vALuE *  •  1 2X *  •  Bf ATHER  INC  CONSTANT  •• EX . 

••VALUE  AT  10  PERCENT* . EX. • INI TI At  V ALUE • */SSX . • C VAPORAY I CK * • 12X. 
••AT  2)  C«/> 

NRITE<  «*733  IOENC  .DC N *. OENX .JEN 25 
7C2  FORMAT (•  *  **0«N3ITY(k3/M2  I*. |6X*KlS.e,3( 13B.PJ8.6U 
•RITE (6.71 |  |v I  SC  *V(SK*V|SA*V IS25 

711  FORMAT!*  *  ,'VISCCSITY  (.'4pAS )  •  ,  1  A  x  ,P  IE.6.3C  IQX.PIS.6)  I 
»RITE(6.73:)  iCL0.4CLK.3CLX 

703  FORMAT (•  •  ,  •  S3LLG IL I T Y  ( G/M3  )  •  •  1 2X *P  IS • A • 2 1 10X .F IS .0 >t 
BRITCI6*  7CA) 

704  FORMAT  <  •  •••V4PCP  ‘•RESSuRf  J  IRE  ACT  CALCULATED  INMRCNOENTL  Y  •//) 
VRITCI6.76S  ) 

70S  FORMAT < *0*  •  *P4CCE S3  CONSTANTS*/) 

•R I TE (6.7Q6IC21 

706  FORMAT  I  *  *  **CVAPC.TAT10N*  *21«.P15.«) 

BR I  TE (6.737 )  CC1.C32.C33 

707  FORMAT!*  *  .  •  C I JPERS ICN *  ,22X  ,F 1 2  •« . 2 ( 10X • P 1 S .2 M 
BAITS  16* 7Ct)CMJ *CM2»CM3 

706  FORMAT  (  •  •  **CMtL4  IP  ICAT  ICN*  ,  10X  .Plf  ,C.  1(  IOX.PlS.au 
BRITEI  6.701  KSl.C32.C43 

704  FORMAT ( *  *  .  •  SRR 5 A  2  I N« • • 2 2 v • P 1 « . I • t I  l C A #F i S • I D 

C 

DO  113  l«i.mtt 
T I M«T I •♦OT  IP 


is: 

124 

12* 


124 
127 
1  2  € 
1  2< 
1  4  C 
14  1 
l«2 
142 
144 
142 
1  46 
14? 
1  4  fi 


14S 

ISC 
12  1 
152 


C 

c 


c 

c 


is: 

1  54 
1  22 
1 54 

15? 
1  26 
lf« 
16C 
14  1 
142 
1  22 


144 
142 
iee 
it? 
i  tt 


1  4S 
1  7C 

171 

172 


1  72 
174 
1  ?f 
1  7 1 
1  77 
1 7e 
l  7S 
lec 


c 

c 


c 

c 


c 

c 

c 


TI*M*T!M/4C  . 
1Ii«*-»T  1MM/60. 
II  MOx TIMm/24. 


CALCULATE  EvA^CPATICR  <^> 

THE  T A  *  TMFTA  ♦  CE 1  *  4 T K * C T I  M  *  1 l  •  “FT  4  }/ VT K 

ftke •  ALOctcxPicip-Af^^cEai-^eiA  «arp.cE2  ♦  u  )/< akp^cs: I 

CVTKE  =-<FTKc-PTK)x  vTK/ (  1  * -FT A  ) 

CFTK  *  MK:  -  FTK 

OVTNE  I^-Tv?N*  (FMAX  -FTM/Cl  .0 -F TK > *£XP < -CT 1M/500 . > 
VTKE-VTKE-  CvTKC 
VTRE»VTN=- CmTNE 


VTEsVTKC*  VTRT 
FCTE=VTE*  t  CC  •  3/V£-*L 
RATKEx-DwTKE/JT  I  M 
RATNEx-Ov  TNE/2T I  M 
RAT ExHl  ATK 2  *R  ATNE 


CtSPERSICN  (C) 

IF  (VlocM.GT.iQOO.)  THEN  CC 
ROKxC.  0 

RDKxCD1mTTK*(*SMS/10.0)m*2.0*CXFI-TTK4<CC2*STK/20.0mCOJmvI3EF/ 
4  10C.C) /Q.CC  1) 

RON  x^  CD  l  *  T  TN*  I  aS^S/  10  •  U  )  mm2  .0  *E  XP  (  -TTN» !  CC2*STK/20  •  0  )/0  *0011 

RONT*RDN*ATN 

RDKT«fiCK*ATK 

DVTKDx-PCKT-OTiM 

CVTNDx-RCKT»CT  IM 

VTKO«VlKD-C VTKC 

vtno* vtnd- c  vtnd 

VT  0*  VTKQAVTND 

PC  T  0» V  TO* 1 CC.O/VSPL 

RATOxRDNTaRCKT 


CALCULATE  MCL3SE  FORMATION 
VI 5«*fcXP ( 2 #2*FM/ (  l-C«l*F»») 
CFm*CM2»!  «S*S*  1 I  4*24 (1 -F«/C*2  )#CT  IM 
Fw  «F •  ♦CF  m 
PC V* 1 00 • • F  • 

Vl  ihM-vl  S»  V I5R 


CALCULATE  SFrEACING  (S) 

CATNS*CJlx { ATN-«0.J3)»EXP< -CS3*T  1N/TTK1  *CT  IM 
CVTNS»OATN$*TTN 

CATKSa  =  °C5  2*  <  I  OO/VIS  CM  »- (A  tK»»  3.  23  )x<TTX-TTK)»*l  .330T  IM-C  VTNS/TTK 


CALCLLATE  Pen  VOLUMES.  AREA.  ETC. 

ATNO«ATN 

ATKCaATK 

ATCCxA  70 

VTK  =  VTK*CVTKE-*-CVTKOACVTKS 
VTNsVTN^OvTREtCVTNJ^CVTNS 
ATNxVTN/TTN 
A  T  K  x  AT  K  ♦  C  A  T  K  3 
TTKaVTK/ATu 


i  a: 

1  «r2 

1  44 
142 
164 
1  €  ? 


ICC 
Its 
ISC 
IS  1 
1  S2 


is: 

I  S  4 

ist 


1  *4 
is? 


I  S£ 
1SS 
2CC 


2  C  1 
2C2 
202 
2  C  4 
205 

iC< 
2  C  ? 


204 
2CS 
21  C 

21  I 
212 
2t : 
214 
212 
214 

217 

214 

2  1  S 
22C 

22  I 
222 

223 

224 

225 
224 


VTC*VTN»VTK 
ATCsATN>  AT* 

CATC*A"C-AT£r 

RAT  jN*<  ATR-ATr.C  )/3T  \t 

RATSKs  (  A T A -1TKC  )  /uT  I  » 

RATjT=(ATJ.at:C)/CTIM 
N I T  *N I T ♦ 1 

C  CALCULATE  n£»  CCMPJSITICNS  OF  SLICK' 

FTNxFMAX 
FT  K*FT  K*DF T  A 

FTx(FTK*VTX*rTR*VTN  )/<  VTK*vTR  ) 

FTNxF  T  N-0  V  lisS-4  IFTN-FTK)/VTN 
TTKCMxTT*/  (  1  .  3-F  M  ) 

C 

C  CALCLLATE  R  6  M  CIL  PROPERTIES 

SCL*SUUU*=xPI-SCLK*FTK  » 

ce r  x  :ENC*(1 .  ♦  :erk  *ftk> 

V  I S«  V  I  SO *P(  V l 'K^F  TK ) 

c  CALCULATE  »:*  CIL  INPUT  (IF  ANT)  ARC  SPILL  POSITION  A?4U  SHAPE 

u 

20R*«SM5*FC3«O T !M 

IF  (V»PL-v:lTo*G.OM  )  220.521.521 
C  5PILL  HAS  S1CPPCC 

52  l  CVCL*  0 .0 

20CxZD0FZCP 
CO  TO  540 

C  SPILL  CONTINUES 

520  CONTINUE 

Z0L«Z0U»Z3P 
s6C  CONTINUE 

Z»x0.5«(-ZCL*(2CL-ZCL*l6.0«ATC/:.l4ie>-*0.5) 

ZL«ATO*4.0/ (2. 14 ie*Z  « ) 

2LE*ZOL«Z«/2 • )*ZCC 
ZTE*2LE-ZL 
C 

C  PP I  NT  TUT  $7  A~G  CF  CIL  SPILL 

C 

IF  (LGUP-NLPJ  52I.52C.520 
520  CONTINUE 

»R  !TC<  4 .5  I  t  >  VSPL 

311  FORMAT (•  l  ».  'VCLU«e  OF  SPILL  (CufilC  METRES)  '.4X.F9.2) 

525  CORTINLE 

•RITE14.5E; )  N IT 

580  FORMAT!*  *  . *  ?  TEA  AT  X  CN  NUM d £ P • .  I  4 X  .  I  I  I ) 

•RITE  <  4.S«J  1  )  UM.TlMM.TlMM.TIMC 

5*1  FORMAT  < '  *  IMc*  .17*. FS.i.iX. •1ECCRCS* »«».FS. 2. 2X,*MniUlES*  .2X.F9 

1  .2  *2X. •HCUPS*  »2X  *F9.2.2X.«CAVS*/  ) 

MR  IT4  <  4. 504  ) 

506  FORMAT ( •  1  .  40  X  »  *  THIN  SLICK*  . 1 6  X •  * THICK  SL I CK • . 18X. 'TOTAL  * . I *PER 
3CE  R  T  •  /  ) 

MRITE!  4.5;?  )  ATN.ATK.ATC. 

507  FORMAT ( •  *  . * AR E A  (SCuAC£  mqT PCS )  •  • 1 5 X • F 1  * . 2 . 2 ( 1 OX.P 1 5. 2 )  ) 

•RItE(4.53C)  VTN.vTK.wTC 

oC9  FORMAT  I  *  *  .  *  V JLuM •  (CU^IC  m? TS * «  )  •  .  J 4X  .  F C. ■  ,  l  3 X .F 1 2 • 2 . 1 3X  .F l 2 .2 ) 

•RITEia.SlS)  ttr.ttk 

5CR  FORMAT  (  '  •.•ThICKScSS  (  M*  T*3E  «  )  •  .  2  7  X  ,F  S  •  4 . 1  t  *  *F«  •  6  ) 

MR  l?£(6.6l  1  )  vTnJ.vTkC .yt-.PCTC 

611  FOR-AT  <•  'MOUNT  I  I 8PL3SCC  (CUE IC  -ETTC « I  • •  12X.F 10 . 5. 2(  1* **F \ C  • 

4SI.5x.F6.: ) 


A-14 


227  4PITE/6.62C tfTN.rTK.FT 

22  6  620  FORMAT/*  '.'FRACTION  2V  APC  -i  A  T  5  2  •  ,  2  2  >  ,F  1  Z  .  6 . 1  SX  .F  1 0  .  o  ,  1  5  X  •  F  1 2  •  5  ) 

225  •RXTE/6.5T5I  *TM1  .VT*C  ,VT  =  .PCTfc 

2  2  C  *55  FC  RMAT (  •  '  •  • A  4C L  N T  E V AP C R A T 2 2  '  .  2 ? X  .  F  1  C  .6  .  I  5 X . F 1 0 .6 . 1 5 A , F  1  0  .6  . 4 X » F 3 

2  .2//  i 

221  »RlTEC6.6lf> 

2J2  eis  format  /  *  '.'process  paTcS  ■»=«  seccss*./) 

232  »RlTt(o,olC)  3 VCL 

2  2  4  ulu  FC«»*AT(*  •  •  *C  I  L  ADO  l  T  I  CN  •  •  57  X  •  r  I  2  •  8  1 

235  »« * TF (6 • 6  l  7 K ATN£  ,HAT<; .RATE 

226  61?  FORMAT/*  *  . 'CVAPCRAT  ICN •  ,  33X.F 1 2  .2  ,  J2X.F 12 . 6  *  1 2X.F 1 2.d) 

227  I  Tt  (  i  .o  t  :  )  R0M,7O«T,R-IT; 

236  ol3  FORMAT/*  •  •  '  C  I  3PER3  ICN • • 24X . P l 2  •  €  .  2 C  12 A *F  1  2 « 6 )  > 

23S  »R1  Tc<  6  .61  «5  )  RAT5N.RA73K .R4TST 

24  C  614  FORMAT/*  '  .  •  SPRE A 31 NG •  , 32X ,F 1 C . A  ,2 I  1 « >  .T 1 C  .  A >// ) 

241  »RITE<o,  334 > 

2  4«  254  FORMAT/*  •  ,  »EV APC«AT ION  CATA») 

242  «R  ITE<  t, J3C  )  VR  AT 

244  J  36  FORMAT/*  ■  . 'EVAPCHAT ION  AIP-CIL  VCLLME  RAT  IC  *  *  16X.F12.3) 

24*  *RXTE( 6 • 3 i 5 >  3TN 

246  J5e  FORMAT/*  *  .  *t;  IL-AATER  INT£RF  AC  l  AL  TENSION*  .24X.F5.2) 

247  »«ITC/6.437>  DEN.CtiNC 

246  437  FORMAT/*  *.»ThICX  SLICK  OIL  D£N 5  1  T V  *  .20X .F 1 C • 3 *6X . • F RESH  »AS»;F|C. 

50 

244  ARlTfc/c.ASe  ISCL.SCLC 

2SC  460  FORMAT/*  *, 'THICK  SLICK  OIL  SCLL B I L  l  T Y •  , 2 1 X ,r \ C. 4 • 6X . • F RE SN  SAS'.F 

610.4) 

261  Afi  X  T  E  <  t  .  4  6  6  )  Vl5.SlSC.VlSF,* 

252  486  FORMAT/*  '.'THICK  SUCK  OIL  PARENT  V  [  SCO  S  ITT*  .  1 5  X  *F  1  C.4.6X,' FRESH 

4  %A  S  *  *F l 0  •  4  »SX  i  *t MuL  S ICN  VlSCCSITT  1S*.FIQ.*> 

252  *R  IT£( 6. 6 2 5  )PC4. visa 

254  625  FORMAT ( *  •.••ATER  CONTENT  CF  TRICK  SLICK  (VOL  PERCENT )•. 1 C X , F 1 0. 4, 

5  10  X*  *VI SC  J  5  I T  V  RATIC».10X,F10.4» 

255  »RITE(e.61C>  TTKEM 

2 56  630  FORMAT/'  », 'THICKNESS  CF  THICK  EMULSION  SL  I  CK  *  ,  l  0X  .  F  1 2 . 5/ / ) 

2*7  4RlTF(t.64C) 

256  C4C  FORMAT / *0  •.  *3°  ILL  LCCATICN  ANC  5F*PS*/1 

2  5S  »R  ITE /  6 .645 >ZL  .*• 

2«C  645  FORMAT/*  •  •  £NCT  H  OF  SPILL  DC  4  NS  I  Nt  •  •  l  2X  •  F  l  2  .  0  •  2X  ,  •  M*  •  5X  ,  •  »  IC  TH  *  . 

47X.F12.0.  .a,*m  •  ) 

261  »RtTE(c*6iC>  CLS.ZTC 

26c  650  FORMAT/*  *. 'DISTANCE  JCSM/ISD*  •  IS*.F 12.0.2X.  *M*  , 5x , • lcaoins  eoce  • . 

•F12.0.PX.**  '  *  *3  X  #  •TRAIL  INC  wDGF  •/  1 

c 

26:  LO  CP  a- 1 

264  *31  CONTINmC 

365  LCCP=LUOP»J 

266  VTKavTK40vCl *3T IM 

267  VSPL»VSPL*OVCL*0T IM 

266  TTK=VTK/AT< 

265  FTKaFTK-C VCL^OT i M«FTK/VTK 

27C  510  CONTINUE 

^  7  1  RETLRN 

2  7  2  ENO 

SCAT  A 
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ILifcC*  I  OQ.OC  CO t*IC 


THICK  'iLlCK  TCI  AL  PERCENT 


90 
9  9# 
999 
-  —  9 
-J  Of  —  <VN 
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90 
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U  J  fOOIMfl 
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2  on  -  - 
K  U  ® 


c  oo 

DN-r 
to  2  IVM39 

u  w  uvrio 

D  —  ♦  9  O  «\!  9  — 

o  3  o#on  —  o 

x  tn  ••*«#• 

<V9OW0« 
2  »-  -  - 
9  V  M 


9  -O 
-OBflfl 
9  O  ♦  tfl  9 
riv»9  w*m 
OfWIUMi 
o#oa,T\i«; 
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a  UIUJ 
*« 
ZU.lt. 

9  9 

o  ui  on 

Off  90 

O  U.  90 

•  90 

«VO  •  • 

-O  9C*0 
3  •  9<0«W 

09  9# 

on  n# 


o  w  vn 

on  90 

o  t,  90 

•  90 

no  •  • 
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at  •  narM 
r*  99 

ON  IV  K 


no*«on 

WOfl^N 

000093 

oroipoi- 


-9  5 

9  OJ  OO 

tua — iu 

a-9  ko 

HIOIIU  4  111 

t*2s£; 

uniat  a  a 

3  —  *  3  *  o 

a  "J — -  >  a 
#3  9*.* 

309—  > 

O-VHJ  2  3 
/l  Ui  W 

—323—3 

1*2-2 
#33303 
*3*2  a# 


3  2 

-  32 

#  -32 
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••22230 
3303 
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—  tt  U 
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u  a  — —  u 
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S  »wu< 
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*  3  2  3  2  U  «  33# 
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•>32*22  *  —3  —  2 
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C>  S^ILL  DCMNn  I  NO  fr  »tlC1l-  246U  .  M 

tlSTAtet  0L»N«lhD  !*!'.£•  *  LfAOfNC.  L'Otit  1*700*  *  THAI  L.  |KC  CCCC 


•  toATPIV  VEH1  LEV5. 

4 

♦  MCNOAY  AUGUST  9.  l  iS  5 


to  Afi  fig  N  STIVER 


MON  I  TOR • VERS ION  3E 
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c 

c 
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c 

c 

c 

c 

c 

c 
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t 

7 


e 

s 

1C 
1 1 
u 
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14 
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2  C 
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22 

24 

2  S 
26 
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C 

C 

c 

c 


c 

c 

c 


c 


c 


tAfiPEK  ST  IV -S  • 


CIL  SPILL  MODEL 
KC*  4  PLEL  OIL 


SET  PA£A*£TEP  VALUES 

SET  COMPUTING  CONDITIONS 
T 1  Mt»  INCREMENT  IN  SECCNOS  -  CT  I  to 
CT  IN* 133. 

TOTAL  TIMc  CC4PUTE0  IN  HCU»S  -  TCCMP 

TC CMP  «  23. C 

MTT*TCC*P«  ’60C/OTXM 

PRINT  CUT  PfifcQOSf.CY  IN  HCURS 

PQPH  *  4.3 

NLP-P0FH4  360C/CT  IM 

MTifl.O 

LOCP*NLP 

SET  SPILL  CCNOITXDNS 

TOTAL  VOLUME  SPILLED  -  VOLTS  IM2) 

VOLTS* 1 00 

DURATION  CP  SPILL  IN  SECCNOS  ( N 1  *  1 <  CNE  INCREMENT)  -  T I MSP 
T 1  MSP* 1 00  •  C 
TI MNP— 7 IMSF/360C. 

T I M*0 • 0 

OVCL*VCLT  S  /  T  1  “4SP 
V! *0 VCL*0  TIM 

INITIAL  ThICK  SLICK  THICKNESS  -  TTK  CM) 

TTK*0.02 

SHEEN  THICK NE  3 S  -  TTN  (M) 

TTN*0. 000)05 
20 L* 0.0 
ZOO*  0 • 0 
VR AT  *0  •  0 

SET  ENVIRCNMF:iTAL  CONDITIONS 
TEMPERATURE  -  TC  ( JEGrtEES  Cl 
TC  *  2. 

TN*TC42T3# 

TP *32  .0»1 .«i’; 

tolNCSPEEO  -  *SKH  IKM/hh) 

•SKM*20 
to$MS*»SKH/3.6 
toSKN* h SM S *  I  .942 
CRIPT  HACTCP  -  FOR 
POP*  0. 0 IS 
VORaFDRtotoSMS 
THETA  *  0 « 


2 1 
2* 
2  C 
21 


3  7 

3  € 
34 

4  C 


42 

42 

44 


4  6 
47 
4  C 
4« 


9  C 
91 


93 

94 


tt 
•7 
99 
S? 
ec 
c  i 


cs 

63 

<4 

69 

66 


6  7 
e« 
64 
7C 


SET  OIL  PKCFCTTXC3 

MOLAP  VOLUME  -  VMQL  CMJ/VCL) 

MAX  PRACTICN  EVAPJHATIO  -  PM** 

ARP  *  0.1  *.  -  10.6  4  <  1.0/203.  -  I  .C/TK  ) 

BNP  ■  30.1  -  AL J G (  1.0  -  C. 00064(14  -  295.)) 

PM  AX  *0.2 
PT  N*0 • J 
PTKaC.C 
P-0.000032 

CIL  RATER  INTCRFACIAL  TENSION  -  ST  CCYNES/CM) 

ST  m  30.0 

«tk«st 

STN-ST 

CIL  OtNSITv  -  CEN J  (KG/M2).  REA  THER  I  NG  CONSTANT  -  OENK 
0EN25  «  ROC. 

CENT-  C.QCC6 

CENC  ■  Og N 2 9  •  <1  -  OCNT * ( TC-25  )  > 

OENK  •  g . | 6 

CIL  VISCOSITY  -  VISO  (CPUI3C).  42ATHER l NG  CCN3TANT  -  VI3K 
VIS2S  a  21.  C 
VI  ST  ■  70GC. 

VlSK  «  5.) 

VlSC-  V 1 329  "EXP( VlST-( | .O/TK  -1.0/299.)) 

CIL  SOLUaiLlTY  -  SOLO  (G/43).  tolATMERI NG  CONSTANT  -  3CLK 
SCLC  *6.5 
SOLK  *  12. C 

CtNX  •  OCNC  4(1.4  OENK  43.1) 

VISA*  Vise  4  2X0  (VISK40.1) 

SOL*  *  SOL C  4  2  XP<  -  SCLK4  0,1) 

9VAPCRAT1CN  CONSTANTS 
CEl*0.C0C674< «SXH)44C. 79 
CE2  4  919. 

CE 3  •  2 J'). 

VTxe-0.0 

VTNgaO.O 

VTE»C.C 

ctSPCRsicN  constants 

CO |40. .Ogci 
CD2-2.0 
CO  3*1.0 
VT*0*0.0 
VT  NO -0.0 
VTO*  0. 0 

6MULSICN  CONSTANTS 

CMl40.es 

CM2  a  2.09-06 

CM3  a  0.7 

PaaO.O 

CPRaO.O 

SPREADING  CONSTANTS 

CSl-t.O 

CSS  *  I  30. 

CS 3*5 . 0 

CSA  *  5  •  0 


A-20 


c 

L  VOLTS  IS  TCTAL  VCLUME  SPILLED  (CueiC  METRES) 

C  VS  PL  IS  V.TLLM2  SPILLED  TC  GIVEN  71  Mg  (CUBIC  METRES  I 

C  TIM*  TIMM,  TIVh.TlMC  ARC  T[M£  (SkCCNCS.  MlkuTiiS.  HOURS*  UAVS1 

C  CT IM  IS  TIME  l NCR TIE NT  (ScCCNOS) 

C  T I  MSP  IS  »P  ILL  CURATIJN  (SECCNCS)  C  MLS  T  CXCEE3  TIM) 

C  OVCL  IS  SPILL  VCLuMc  INCH£*CNT 

C  PM A A  IS  MAXIMUM  FMCTI-N  CF  The  ClL  HHICH  CAN  EVAPORATE 

C  ttSKH  •  ASMS  AHJ  A  SAN  **".6  4  INC  S°EE3S  IN  KM/H*  1/S  AND  KNOTS 

C  FOR  IS  »(n:  JUFT  FACTLR 

C  TTN*  T  TK  »  TT0  ARC  ADS  Th  I N  •  THICK  ANC  AVCPAGE  SPILL  THICKNESSES 

C  ATN  *  ATK*  ARE  THIN,  THICK*  A  F  3  TCTAL  SPILL  AREAS 

C  ( SOLAR L  Me TrfCS  I 

C  VTN.VTK,  V  TC  ARE  7H I N,  THICK  AND  TCTAL  SPILL  VOLUMES 

c  <cuaic  McrpES) 

c  FTN.  F  TK  AND  F’  A+Z  ThIN,  THICK  AND  TCTAL  FKACTICNS  EVAPCR»T  =  C 

C  ST,  STK.  AND  STK  AP£  INITIAL,  AND  THIN  ANC  THICK  SLICK  INTERFACIAL 

C  TENSIONS 

C  CENG.  V I  SO  AN  J  SOLO  AK  £  INITIAL  CENSlTV*  VISCOSITY*  AM)  SCL'JOtLlTY 

C  CF  THE  CIL 

C  DENT,  VIST,  V"*RT  AH  E  CCNSTANiS  RELATING  DENSITY*  VISCOSITY  ANC 

C  VAFCR  PRESSLFE  TO  TUMPEPATURE 

C  OEN2S*  V S  32 S •  VPRJ5  APE  FAtSH  CIL  CENSITY*  VISCOSIT <  AND 

C  VAPCR  PREiSLPC  AT  25  DEGREES  C 

C  CE  NX •  VISA,  S  CL  X  ,  VPRX  APE  VALUES  CF  DCNSITV*  VISCOSITY  AND 

C  VAPOR  PRESSURE  AT  ICX  *EATHERINC 

C  INITIAL  RAT  INDICATES  FATES  CF  PROCESSES 

C  ZL.  Zw  ARE  LENGTH  ANC  «IDTH  CF  SLICK 

C  2DL  IS  DRIFT  LENGTH  CF THE  SPILL  SOURCE  POINT 

C  200  IS  DRIFT  LENGTH  FROM  LAST  SPILL  PCINT 

C  TE*P  IS  EXPOSURE  TIME  (SECONDS) 

C 

C  INITIAL  0  MEANS  DELTA  CR  DIFFERENCE 

C  FINAL  E  INDICATES  EVAPORATION 

C  FINAL  C  IMCICATES  dispersion 

C  FINAL  M  INDICATES  EMULSICN 

C  FINAL  S  INDICATES  SPREADING 

C  FINAL  »  INDICATES  DRIFT 

C 

C  SHUT  orF  FACTJRS 

7 l  SOFFE« I  .0 

75  SOFFQ*  )  • 0 

7  z  SQFFMa 1 *0 

7o  SOFFS* 1,0 

7  S  SOF  N*  I  .0 

C 

C  SHUT  OFF  PROCESSES 

7  t  FMAX=F MAX* SCFFE 

77  C  E  l  *  S  OF  F  E • C  6 1 

7  €  CD  1  *  SOF  F  O  ^  C  C  l 

7S  CN2»S0FF4*CM2 

El  CSI«CS1 *SCFFS 

€1  CS2«CS<!  •  3GPF  5 

e<  F0R»F0h*5CFFM 

C 

C  INITIAL  CONDITIONS 

c 

c  calculation  cf  initial  disposition  cf  cilslick 

B2  VSPL*VI 

ex  ATK*VI /(T T KtCSAXTTN  ) 


fe* 

et 

£7 

Ef- 


SC 
9  l 

5* 

s; 

94 


St 

St 

712 

S  7 

St 

71  J 

sv 

1  Ct 

71  7 

1  C  i 

1  C  5 

715 

C3 

IGA 

JSC 

:  os' 

ice 

351 

l  C  7 
ICE 

710 

ICS 
*  1C 

700 

1 1 1 
l  12 

701 

1  12 

1  1  4 

702 

1  IS 

1  It 

7  1  1 

1  1  7 

1  1C 

0  3 

US 

12C 

704 

12  1 

125 

70S 

122 

124 

7C6 

12! 

1 2 1 

707 

127 

I2C 

706 

129 

1  J  C 

7-9 

12  1 
132 

C 

ATN«ATK*C  i4 
atc*atk.a~n 
VTKsATm  Tin 
V  TN*  ATN  •  T  Tn 

calculate  initial  oil  PROPERTIES 

CE  NS  IT  r 

CEN  -  CENC ♦  (  I  ♦  ♦  JENK • FT  K  ) 

vi scasi rr 

VIS  «  Vise  •  EXP< VI SKXFTK  ) 

VI ScM*VJ  S 
SOLUBILITY 

SOL  =  SOLO • E X  P ( - SCL  K*F  TK ) 

•R ITE(C,5CS>  VCL73 

S  FORMAT ( »| *, 'TOTAL  VOLUME  TO  QE  S P I LLEO* ' . F 9 . 2 , 2X * • CUD I C  METRES*/// 
l  /✓) 

1R1TE(C.715)  TIMJP.TIMHP 

FORMAT (  ■  •  ,  *OURAT  ION  OF  SB  ILL  ( S  )  '  .12X*F9,  I  •  SX  * ' OR •  *Fll,3*2X**M*  ) 

XRlTt (b ,71  2 )  V  I 

FORMAT (  1  *  ,  *  j  »  ILL  VOLUME  INCREMENTS  (M3  )  •  ,*x  ,F9. I ) 

vRITE(c#T;7»  TTK.fTN.CSA 

FORMAT  (  1  ••'lUTtu.  SPILL  THICKNESSES  (  M  >  •  •  3X  *F  I  |  .6*  2X«  •  TH  ICK 
SL  ICK* ,sx, F  1  1.6.0X.  • THIN  SLICK*  *«x  ,F  11.6, 2X.*ARZA  FACTOR*//) 

■  RUE  <6.71  «  ) 

FORMAT! *C»  ,  *FN VIE  DNMcNTAL  CCND ITICNS*/) 

•R  ITE( 6*  33C )  TC  *TF 

FORMAT! *0*  .  'TEMPERATURE  <  DEGREES  C )  *  *  1 0 X . F9 • 1  * SX * *0R • * F 9  .  I • 2 X * 

‘  *(  0EGREC5  F  )  •  ) 

»R  ITE(6*35  1  )wSKH*rfSMS*MGKN 

FORMAT!*  •  *  * •! NC  SPECS  ( KM/H  )  *,  1  EX  •  F 9 . 1 • SX •  *0R *• F9 • I , 2X , *M/3  •* SX • 

1  *0R*  *F9* I • IX*  *  KNCT3 *  > 

»RXTE(6,7lC>  F  CF 

FORMAT!*  *, 'DRIFT  F AC TCR *  .2 6 X  ,  F £  .  J / / 1 
•RITE  (6. 7  30 

FORMAT! *C ». *CIL  PROPER T I ES  '/ I 
•RITE! 6.731) 

FORMAT (• c 1  .* INI T IAL  VALuE  •*  1 2X *•  »EA TMEP  INJ  CONST ANT *. 6X * 

•VALUE  AT  10  PERCENT*  *oX.* INITIAL  VALUE  *  */86X« ‘EVAPQHAT ICN* •  I2X* 
•AT  25  C*/  I 

•R!T2(6«732  )OCNC  ,  Jc NK.S0NX.JEN25 

FORMAT! •  ••*J&N3ITr(XG/w3>>*16x*FlS*6*3(lOX*FlS*6>) 

•RITE (6 *7i 1  )VlSC*VtSK*vISX,VlS^S 

FORMAT!*  •  ,  *  VI SCC31 TV  ! MPAS )'  ,14**7  1 3  *6  •  J  (  l  OX *F l S. o  ) ) 

•RtTE(6,73:i  iCLO  •  30LK  *  SOL  X 

FORMAT ( *  •  ,  *5 JLUCIL JTV  ( G/ »3 ) *  * ) 2X ,F IS .6 , 2 ( I  OX , F 1 5.6 ) ) 

•RITE! 6.75A  ) 

FORMAT!*  *.* VAPQR  PRESSURE  3  AR  2  NOT  CALCULATED  I  NOE  PENOEMLY  •  ✓/) 
•RITE(6.703> 

FORMAT!  *  3  '  .  *P'-1CLt  33  CJNSTANT$*/I 
•RITCI6.7 )C“ 1 

FORMAT!*  •  ,  *lVAPCSAT ICN*  .ZlX.FU.fl) 

•N ITt!6.727)  CC1.CD2.CJ3 

FORMAT  (  •  •  ,  ‘CISPS^S  I CN  *  *  22X  ,  F  l  •  ,  fl  .  2  <  10X.Fl?  .CM 

•RITE  I  6. 7CE )CM1 .CM2 .CMC 

FORMAT  I  *  *  *  •«  4UL5IF I CAT ICN*  *lEx*FIS,e*2(lCX.F12.C)) 

•R1TEI6*  709)Cit  .C32.CSJ 

FORMAT!  •  •  #•$.*  1CAJINC*  *2  3  X  ,F  |  •  ,  C  ,  ?  (  |  CX  .F  l  S  .  a  )  ) 

CO  510  L*l  .NITT 
TJHmTImHOT  IP 


A- 21 


71 »»«TXM/«C. 
T|M*«Tl«*/eC 


132 

134 

138 


128 

137 

138 
Uf 
14C 
141 
145 
143 
1  44 
148 
1  4t 
1  47 
14k 


144 
ISC 
18  1 
152 


C 

C 


c 

c 


152 

i  54 
16  8 
i  te 
157 
156 
18  4 
16C 
161 


1  6  4 

165 

166 
167 
1  66 


164 
I  7C 
17  1 
17k 


1  72 
1  74 
1  7£ 
1  76 
1  77 
176 
17V 
1  6  C 


C 

c 


c 

c 
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CALCULATE  »4A3C»ATX0N  (51 

7HtTA  a  THETA  ♦  C  C l  •  AT* ♦  CT  X  *»  •  (  1  .-5  TK  1/ VTK 

5T*t  a  ALQC (ExP( 0 fl a- C53 > * T ► 5 7  A  •ARP4CE3  ♦  1  •  )/ (  ARP4C63  ) 

CV7K.E  •-{P7KE-5TKJ4  VTK/ <  1  ,-5T*  1 
05TK  a  5T*C  -  5TK 
VfiATa  THE  7  A 

CvTNE  ■  -  >TN*(5MAX  -  5T  N  )  ✓  (  1  •  C  -  6  7  K)  •  IX»»  (  -  CT  I  M/80  0  •  ) 

VTXE»V7KC-CVTKS 
V7NC»VTNt-CVTNE 
VTiavTKC^ vTNE 
PC7EaVTE«l CC. O/VfPL 
XATKEa-DVTKS/3TXM 
PA7nE«-OvTN£/07:n 
AA  7EaRATKE+RATNE 

018PERSICS  (C) 

15  (V13CM,GT. 1 CCC.)  THEN  00 

RDK-O.C 

ELSE  DC 

RDK-CQ 1 ■TTX4(KSMS/l0.0)442.04C*O(-T7X4(CD2*STK/i0.O^CO34Vl 856/ 

>  100,0) /O.OC  1  ) 

END  15 

RON  ■  CDI-TTN4(nSMS/10.C)m4?!.C*E>P(-T7N4(CC2*3TK/20,0)/0»OC1 1 

RDNTaRON*  4  7N 
RDK7«RCK-ATK 
CVTKDa-RCXT40Tll# 

CVtND— R0N74DT  IN 

VTKD*VTKO-CV?KC 

VTN0*V7N0-CVTNC 

VTD-VTXO+VTNO 

PCTC-VTD- l CC.3/V5PL 

RATD-RONT^RCKT 

CALCULATE  *CuSSC  5GR  X  AT  I  CN 
VI £A-EXP<  2 .5"5»/< 1-C5145*) ) 

CF»«CM2»(  *  a*S*  1  )*42Ml-5*/CM2)*CTlN 

P*-F  »♦ uF  • 

PC-*  100.4F  . 

Vl SEM-V13- V 1S» 

CALCULATE  SPREADING  (S) 

CATNS=»t  si  •  (  A7Nt*0«33  >»£XP  (  -  CS3  ■  T  TN/ TTK  )  *0T  1M 

CVTNS-0A7N<-TTN 

CV  7k  S*  ~D  V  T  N  * 

CATkS  *  C5?- t 1 00/ v I  SEN) 4 ( ATK--0 .23  )  • (TTK-TTN) •• 1 •  3J*3T Z  M-D VTNS/TTK 


CALCULATE  N£»  VOLUMES.  AREA,  ETC. 
ATNO-A  TN 
AT  KG*  AT  K 
ATCCsATO 

vTKivTK-OvTKE-CVTO-CViKS 
VTNsvTN*nv7NC4CvTN0*CV7N3 
AT  N= VTN/ T  7N 
AT  K» AtK«OA TKS 
7TK*VTk/AT  •« 


18  1 
112 
112 
1*4 
118 
1«< 
117 


1  •« 

188 
1  8  C 
18  1 
18# 


183 

1  84 

188 


1  91 
IV? 

184 

188 

ICC 

20  1 
261 
d  02 
2C4 
208 
2ce 

207 


2ce 
*ct 
21  t 

21  I 
21  i 
2  1  2 
214 
211 
516 


2  1  7 
216 


218 
2  2  C 
221 
22? 
12  J 
224 
22  8 
kit 


C 

c 


c 

c 


c 

c 

c 


c 


c 

c 

c 


VTC-VTN* VTK 

AT  C-ATN+  AT X 

OATC«A7r.-A7CC 

RAT  SN- I  AT N -  AT  <Q ) /CT 1 M 

AATS»K«(  ATK-ATKO/OT  IV 

«AT17a( ATO-ATJC) /0T1 5 

MTaNlTM 


CALCULATE  NE4  COMPOSITIONS  05  ILlCMl 

5TN-5MAX 

5TX-F  7  n ♦ Or  7  N 

57 «(fTN-VTX«PTN*vTN)/( VT*4 VTN ) 
5TN«57N-OvlN5*(67N-5TX)/V7N 
TTKEM-TTK/ (  1,0-5*) 


CALCULATE  *C*  CIL  PROPERTIED 
ILL- OUL Z >P(-SCLX45  7K  ) 

CEN  •  CENC  •  (  l  .  ♦  30*  •PTX) 

V 1  8-  v  1  SC-c  JT(vI3K«5T*> 

CALCULATE  Ne-  CIL  K-PwT  <  15  ANY)  AND  SPILL  PC31T10N  A  40  SHAPE 

2DB-*0vS-5CR^ JT JM 
15  (vSPc-V  JL^fC.CCl  )  825.821,53  1 
•PILL  MAI  S  7  L'PPf  C 
111  CVOL-C.O 

20C-20C42CR 
CC  T r  560 
SPILL  CONTINUES 
510  CONTIME 

ZDL-20L+  2U* 

560  CONTINUE 

Z»-3.S»(-2CL*( 2CL*2CL«16 «0«ATC/3 . 14  U ) ••0.5) 

2L-ATU-4. 0/(3,141 6* 24) 

2L  E«20L42«/2 *0 4ZC0 
27E-2I.E-2L 

PRINT  CUT  STATE  C 5  C  IL  SOJLL 


If  (LOCP-NLP)  521.530,530 
330  CONTINUE 

*»1 Tfct 6 .51 1>  VSFL 

sit  FORMAT (• 1 •, 'VJLUMC  OF  5P  ILL  (CufIC  XETRE S )  '  ,4X ,5 9,3 ) 

535  CCNTINoC 

■KITF<c.3J:)  nit 

510  50NJ*  AT  {  *  *  •  *  I  T  EM  A  T  I  ON  NO  MU0  P*.  144.1  11) 

•  PITt  (c.fitll  )  T  IN  ,T!  M  ,T  INH,-  INC 

661  FORMAT  (  •  •  .•"'INC*  ,17X,6  9,?.  2X.  •  «:CCK3S’  .2  k  ,5  7. 2  .2X  .  •  NI  NU  T2  S  *  .2X.  59 

1  .1  .iX, •MCJPO' «*X  ,54,2.2k,  ' OATS  1 /  ) 

•r i rt ( t . s:e ) 

aOS  5C  AT  (  *  •.-CX.'TmIn  SL  1  CK  •  .  I  6  X  ,  •  Tm  I CK  SL  1  CK  •  ,  |  0X  .  •  T  OT  AL  1  •  1  3  X  •  •  *»£R 
3 CENT* /) 

4RMF(£,sC7)  ATN.4TX.ATC 

aCT  5CR/ATC*  •  ,  • AT  E  A  (SOOAPG  “27 s« 1  )  •  ,  1  8 X  ,F 1 5. 2 , 2<  iOX ,5 1 5. 2 )  ) 
»R!T0(6«SJt)  vTN.vTK.vTO 

50rt  5(.RxaT(*  *.*VLLUMr  (  c  v/tf  X  C  V  •  7  P  |  )  •  .  2  4  k  ,5  4 . 5  ♦  I  OX  .  5 1  2 . 2  •  1  3*  •  5  l  2  *2 ) 

*5 ITL (6. IDS  )  ’TN.TTK 

008  5CM'AT<»  *  .  *T<:CXNCS«  (MtTBHS)  •  .27k  ,X9,6, 16X.59.6) 

4N  175(6.61  tl  VTNC.VTXJ  .VTO.PC T3 

oil  FORMAT  (•  *  .  •  ANCLNT  CISPF.PkCC  (CufilC  vgTFCS  »  •  .  12X.5  10.3.  J(  ISA.5X  C. 

48) .5k. 6 1.2  ) 


A-22 


22?  *RIT6{6.62C  )FTN,FTK,FT 

2*e  620  FORMAT! •  *. 'FRACTION  E V APCR *T * C  •  ,23 * «F 1 3  .6 . 1 5X . F 1 0 . 6 . 1 5 X , F l 2  .a  I 

22S  »RITE!6.S3* >  VTNC  tVTKC  .  VTf=  #  i»C  Tr 

23C  582  FORMAT!*  *.'A»CwNT  6 V *p C R AT C C •  .  2 7* • F I 0  .6 .  1 3X . F 1 0 .6 . 1 5X , F  l  3 .6 • 4X. F3 

2  .2//J 

22 1  *R I TE  (  6  » A  1 1  ) 

221  fe  1 5  FORMAT!  •  •.  'P^CCCSS  5AT«?g  pcr  S5CCNC*/) 

222  4HlTE(c.616)  jvjl 

234  ot6  FORMAT!*  '  A  CO  I T  I CN  •  ,  57  X  ,  F  1  2 . 8  ) 

22  5  *RITfc<  c.61  l  IRATNE  .aATK£  .  A  T  = 

236  PIT  FORMAT!*  •  .'EVAPulAT  ION*  .3.3A.r  12,6  .  I  2 X  ,F  1  2  .  «  .  I  3  X  ,F  l  2  •  6 ) 

23?  •RITLC6.MS)  H CNT  . 3C * T  . O AT c 

232  alrt  FORMAT!*  •  .  •  O 1 $PF AS  I  UN •  .  3* A  .F  l  2  •  C  •  2  (  13 A *F l 2 . 9  )  ) 

23R  4R  !  TE  !  6  1 6  1  *>  l  3  AT  i.l.  ft  ATS*  *  RAT  3T 

2  AC  CIS  FO  KM AT (  *  •  •  •  $ TRg A 3 t NG •  •  3 3X  .  F  l  0  .  4  ,2 (1  «X . F l 0  .4 )// ) 

241  *R ITE ( o«324) 

242  354  FORMAT  { *  •  ,  * E V APCRAT  ICN  CATA') 

242  »RITE!6.32t)  Vfl  A  T 

244  356  FORMAT  1  *  *  .  '  E  V  4PC-1AT  ION  AIR-CIL  VOLL«E  RATIO*  »18X*F12»3) 

245  4«ITe<o.355>  3TN 

24t  J5d  FORMAT  <  *  »**ClL-4ATcS  INT£RFACIAL  T  E  N5 1  CN  *  ,24*  .F  5.2  i 

247  4RlTt(c.437)  CEN.OENC 

2  4  €  467  FORMAT!*  *. 'THICK  SLICK  CIL  06 N 5  it Y  •  •  20X • F  10 . 3 • oX • • FRESH  MAS*. Ft  C. 

50) 

24S  »RIT6 (6  .465 >£CL.SCLC 

2  2  C  48d  FOWMATI*  •••THICK  SLICK  CIL  SCLLc  I  l  JTV*  , 2 » X  .F t 0 .4 »6X . • FR63F  6AS*»F 

6  10  .4) 

221  •RtTC(6»4d6)  w  l  S  •  V  I S C  ,  V  I SE N 

222  466  FORMAT!  •  *. ’THICK  SLICK  OIL  PARENT  V  I  SCC3 IT V •  » 1 5X »r i 0 .4 .6X  .  • F RES F 

4»AS*  .F10 .4  .5X,  'EMUL3 ICK  VISCOSITY  iS*.rl0.4) 

252  4R1T6I6. 625  )PC4. VISR 

254  025  FORMAT!*  *.*MATER  CONTENT  Or  THICK  SLICK  ( VCL  PERCENT )••  1 0 X  ,  F 1 0* 4 . 

510*.*VISCCSITV  R4TlC*.l0X,Ft0.4) 

255  •RITEle.tJC)  TTKFM 

2ht  o30  FORMAT  (•  *  .  •  Th  1  CK.NES 5  CF  THICK  EMULSION  SL I CK  '  •  1 0X  .F 12*  5//  ) 

22?  *R  I  T  £  !  6-  «  6  4  C  ) 

256  640  FORMAT  (  *0  *. 'S  MLL  LOCATION  ANC  $»-ARE*/> 

25R  »R  ITt  1<  .04  2  )  2L  ,24 

2<C  c  4  5  FORMAT!  •  *. 'LENGTH  CP  SPILL  3 C «N 4  I NC •  » l 2X , F 1 2 . 3 • 2X . * M* . S A .  * *  I CTM * . 

4  7X.F 12  .0.2*  .  *  1  •  ) 

261  4RITEI6.62C)  ZLC.2T6 

262  o53  FORMAT!*  •.•DISTANCE  JC4 N» I  NO •  •  1 R X  .  f 1 2 • 0 *2X  •  •  A •  . 5X • • LEAD  I NG  EDGE*. 

*Fl2*0»2X«'M*.3Xf  'TRAILING  cUGE  •/  1 
C 

26:  LOCPc-l 

264  531  LCnTJnlc 

265  LCCP*LC0P>1 

2  6  t  VTK  =  VTK »0VCL*3T  IM 

26?  VSPL ■ V CPL ♦ 2 V CL "CT  Im 

266  |T K« VTK/A T  K 

2 6<>  FTK*FTK-3vrL»0T  IM4FTK/VTK 

2  7  C  SiO  CONTINUE 

271  RETURN 
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